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ABSTRACT 


TIiIb  report  consists  of  the  derivation  of  approximate  analytic  solutions 
to  the  interior  ballistic  equations  of  cartridge-actuated  devices.  A  unique 
feature  of  the  system  is  that  linear  burning  rates  need  not  be  assuned.  The 
solutions  are  presented  in  graphical  Tom,  thereby  pennitting  the  rapid  solu¬ 
tion  to  a  variety  of  ballistic  problemsr  Several  sample  problems  are  included. 

A  complete  list  of  the  types  of  problems  discussed  is  as  follows: 

A.  Detemination  of  maximum  pressure  and  ejection  velocity  for  given  pro¬ 
pellant  loading. 

1.  Linear  form  function 

2.  Quadratic  form  function 

3.  Horizontal  and  nonhorizontal  firing 

4.  With  and  without  covolxjme  correction 

B.  Detemination  of  grain  design  parameters  to  yield  given  performance 
requirements. 

1.  Given  maximum  pressure  and  ejection  velocity 

2.  Given  maximum  pressure 

3.  Given  ejection  velocity 

C.  Detemination  of  ballistic  effects  of  small  changes  in: 

1.  Mass  accelerated 

2.  Impetus 

3.  Charge  weight 

4.  Web  size 

5.  Burning-rate  coefficient 

6.  Chamber  volime 

7.  Stroke  length 

8.  Cross-sectional  area  of  tube 

D.  Design  of  a  scale  model  to  reproduce  balljntics  of  a  larger  cartridge- 
actuated  device. 
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INTTRODUCTION 

In  experimental  work  with  catapults  and  other  cartridgeractuated  devices, 
it  is  frequently  desirable  to  be  able  to  predict  the  maximum  pressure  and 
ejection  velocity  for  a  given  loading  with  a  minimum  of  computational  effort. 
The  interior  ballistic  performance  of  a  cartridge-actuated  device  is  governed 
by  a  system  of  non-linear  differential  equations  which  can  only  be  solved  by 
tedious,  time-consuming  numerical  methods.  Indeed,  a  single  solution  by  such 
methods  often  requires  as  much  as  eight  or  ten  hours  using  a  desk  calciilator. 
Although  many  short  methods  have  been  developed  for  the  determination  of  maxi¬ 
mum  pressure  and  ejection  velocity,  (b-g)  *  the  assunptions  upon  which  they  are 
based,  e.  g. ,  linear  burning  rates,  neglected  energy  terms,  etc. ,  quite  often 
lead  to  errors  of  considerable  magnitude  in  cartridge-actuated  device  applica¬ 
tions.  .  • 

It  is  also  desirable  to  be  able  to  solve  the  ccmpanion  problem,  i.  e. , 
the  charge  weight  and  grain  dimensions  required  to  proiuce  a  given  maxiimxn 
pressure  and  ejection  velocity.  As  far  as  is  known,  until  this  time  no  npn- 
errpirical  method  has  been  developed  for  solving  this  problem  directly,  it 
being  necessary  to  employ  a  nonconverging  trial-and-error  process. 

The  purpose  of  this  report  is  to  provide  a  simple  but  accurate  solution 
to  both  types  of  problems.  A  typical  calculation  using  the  methods  developed 
here  can  be  completed  in  about  15  minutes  and  requires  nothing  more  than 
training  in  high-school  mathematics.  Resiilts  appear  to  agree  to  within  5  or 
6  per  cent  of  those  obtained  through  numerical  integration,  although  an 
exhaustive  comparison  of  results  has  not  been  made.  Simplifying  assumptions 
are  introduced  which  permit  an  analytic  solution  to  the  ballistic  equations. 
The  results  are  presented  as  families  of  curves  relating  certain  dimensionless 
variables.  Numerical  examples  are  then  presented  to  illustrate  the  use  of  the 
curves  in  solving  the  types  of  problems. 


^Letters  in  parentheses  indicate  references  at  end  of  report. 
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The  nomenclature  used  berQ  has  been  made  to  eorrespand  as ‘closely  as  pos¬ 
sible  to  that  enf>loyed  .in  a  previous  Atlantic  Research  Corporation  report  pub¬ 
lished  in  1958  (i) Particular  attention  sViould  be  given  to*  the  units 
assigned  to  the  parameters.  These  are  included  along  with  the  nomenclature  in 
Appendix  G.  '  ‘ . 

For  brevity,  a  cartridge  actuated  device  is  hereinafter  referred  to  as  a 

CAD. 


I.  HJNDAMEMTAL  EQUATICIJS 

The  equations  presented  in  this  section  are  for  the  most  part  derived  in 
the  report  (i)  cited  above.  They  are  included  here  to  emphasize  changes  in 
notation  and  to  make  this  a  self-contained  report.  The  basic  equations  are 
divided  into  two  g^ups,  those  applicable  during  propellant  burning  and  those 
applicable  after  the  propellant  has  been  consumed  (after  "burnt V) ,  ■ 

A.  During  Burning 

1.  Tte* Equation  of  Motion  of  .the  load. equates  the  net  force  on  the  pis¬ 
ton  to  the  mass  in  motion  times  the  acceleration,  as  follows: 

•  ■  ■  (1-1) 

where  is  the  pre'ssufe  on  the  base  of  the  piston,  A  is  the  cross-sectional 

area  of  the  piston,  Ff  is  the  total  resisting  force  due  to  friction,  u  is.  the 

gross  accelerated  ma^s,  g  is  the  acceleration  due  to  gravity,  --  is  the  fac- 

.  ic 

.  *  •  * 

tor  for  converting  pounds  of  mass  to  slugs,  -9  is -the’  angle  "Of  launch  measured  . 

from  the  horizontal,  x  is  the  distance  traveled  by  the  piston,  and  t  is  the 
time  measured  fron  the  beginning  of  motion  of  the  load. 

2.  The  EguatiOT^^f^Sta^  is  the  relation  between  the  pressiare,  density,  • 
and  tenperature  of  the  propellant  gases.  The  Abel  equation  which  applies  with 
sufficient  acciiracy  for  most  CAD's  may  be  written, 

r  12  Ax  +  4/^  ,  1  1  ,  ^ 

P  1  - ~C - iv  -  ^)  j  ^  12  NRT  .  •  (1.2) 
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irtiere  p  is  the  space-average  gas  pressure,  is  the  initial  free  volume 
available  to  the  gases,  Q  is  the  mass  of  propellant  burned  at  time  t,  77  is 
the  covolume  factor,  p  is  the  density  of  the  solid  propellant,  9  is  the  num¬ 
ber  of  moles  of  gas  per  unit  mass,  R  is  the  universal  gas  constant,  and  f  is 
the  space-average  gas  temperature.  The  factor  12  is  introduced  to  correct 
for  the  inconsistent  units  enployed  here  as  in  most  interior  ballistic  sys¬ 
tems. 

3.  The  ]ft>ergy  Balance  Riuation,  assumicg  no  gas  leakage,  states  that 
the  total  energy  lost  by  the  propellant  gases  is  equal  to  the  kinetic  energy 
of  the  load  in  motion,  plus  the  potential  energy  of  the  load,  plus  the 
kinetic  energy  of  the  propellant  and  propellant  gases,  plus  the  energy  lost 
through  heat  transfer  through  the  CAD  walls,  plus  the  energy  dissipated  in 
overcoming  friction.  Thus, 

-  rt  =  (^)  '  -  fy-o  (1.3) 

where  is  the  average  constant- voltnne  specific  heat  of  the  propellant 
gases  over  the  range  of  teirperatures  in  the  CAD,  Ty  is  the  constant  volume 
flame  temperature  of  the  gases,  Bp  is  the  kinetic  energy  of  the  propellant 
and  propellant  gases,  and  Ef,  is  the  energy  lost  through  heat  transfer  through 
the  launcher  wall. 

The  Burning  Rate  Equation  of  the  propellant  expresses  the  propel¬ 
lant  burning  rate  as  a  function  of  the  space-average  pressure.  The  follow¬ 
ing  equation  fits  the  experimental  data  of  most  propellants. 


r  =  Bp”  (1.4) 

di stone© 

where  r  is  the  propellant  burning  rate  with  dimensions - ,  n  is  the 


burning- rate  exponent  with  a  usual  range  of  0.2  <  n  <  1. 2,  and  B  is  a  con- 
stant  of  proportionality. 

•  5.  The  final  relation  required  is  the  Form  Jtoction  of  the  propellant 
grain.  When  a  solid  propellant  bxims,  all  the  uninhibited  surface  recedes 
at  the  same  rate  assviming  simultaneous  ignition  of  all  surfaces.  Thus,  there 
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is  a  single  geometric  relation  called  the  form  function  which  relates  the 
mass  of  the  propellant  burned  at  any  time  to  the  distance  burned- through  the 
grain  at  that  time.  Most  propellant  grains  presently  in  use  have  a  cubic 
form  function  of  the  type 

C  =  A,i  +  A2  +  A3  (1.5) 

where  Aj,  A21  Ag  are  constants  depending  upon  the  geometry  of  the  grain,  and 
L  is  the  distance  burned  through  the  grain.  (See  Appendix  A  for  a  discussion 
of  the  form  functions  of  some  of  the  more  CCTimon  grain  configurations.) 

B.  After  Burnt 


Let  2  be  any  one  of  the  variables.  The  value  of  2  at  grain  burnout  will 
be  denoted  by  will  be  referred  to  as  the  "web",  *  and  as  the  charge 
weight.  After  burnout,  the  ballistics  are  governed  by  equations  (l.  l) ,  (l.S), 
and  (1.3)  where  C  is  replaced  by  the  constant  C^. 

II.  DEVELOFMFOT'  OF  THE  FQUATIONS  OF  MOTION  AND  ENERGY  BALANCE 


A.  Kinetic  Energy  of  Propellant  and  Propellant  Gases 

The  motion  of  the  propellant  and  propellant  gases  is  described  by  the 
usual  partial  differential  equations  of  fluid  flow.  The  solution,  however,  is 
entirely  too  cumbersome  to  be  included  in  any  interior  ballistic  system.  An 
approximate  solution  developed  by  Kent^  is  used  in  most  systems. 


*  Some  interior  ballisticians  define  the  web  as  the  minimim  distance  between 
the  initial  sui'faces  of  the  grain.  Using  this  definition,  the  web  is 
usually  equal  to  2Lj,.  Care  should  be  taken  so  as  not  to  confuse  the  two 
definitions. 

^See  reference  (b),  pp.  141-147  for  the  derivation. 
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According:  to  this  solution,  the  spac^average  pressure  is  related  to 
tbe  pressure  on  the  base  of.  the  piston  by  the  egression 


(2.1) 


where  cr  is  a  constant  depending  the  ratio  C^/w,  with  a  \isual  value  of 
about  3.  The  kinetic  energy  of  the  propellant  and  propellant  gases  ia  related 
to  that  of  the  accelerated  mass  by 


(2.2) 


The  value  of  the  ratio  cy<ru  for  most  CAD'a  is  of  the  order  of  ^nce, 

little  inaccviracy  is  introduced  by  setting  .=  i>  in.  equation  (l.  1) ,  and 
Ep  =  0  in  equation  (1.3). 

B.  gr jetton  Biergy  Losaea 

The  friction  force,  Fy,  ia  assumed  to  be  a  constant  percentage  of  the 
pressure  on  the  base  of  the  piston,  Pg.  (See  reference  (b),  pp.  8-9  for  Justi 
fication  of  the  assiinptlon. )  Let  the  constant  of  proportionality  be 

K  A 

,_1 —  .  then,  since  P-  .=  p,  equation  (l. 1)  becomes 
1+f,  ’  S'''* 


PA 


(2.3) 


The  assunption  of  proportionality  also  permits  the  evaluation  of  the 
integiral  J  F^dx  in  equation  (1.3).  Let  v  .=  g,  then 


/  '  '■>7  *  2  »»  «•  (2.4) 

0  .  .  0  .  *  *'  • 

Setting  (l  t  jr^)'  w  =  Wj,  an  "effective  mass"  being  accelerated,  equations 
(1.1)  and  (1.3)  slnpllfy  to 


(2.6) 
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and 


CyC  (Ty  -  T)  =  (v*  +  2  Bin  0  )  +  £/,  . 


(2.6) 


The  value  of  is  usually  obtained  eirpirically  from  eject  ion- velocity 
data,  and  an  integration  of  experimental  pressure- time  curves.  It  follows 
from  equation  (2.3)  that 


J  pdt  -  (v^  +  g  t„  Bin  0) 

a 


(2. 7) 


where  subscript  m,  here  and  elsewhere,  denotes  values  at  the  ccirpletion  of  the 
piston  stroke.  Rearranging  equation  (2.7),  the  value  of  A'j  is  foimd  to  be 


Ag, 


,/■> 


pdt 


K,  = 


*  +■>■  t'"Tin"  0)  ■  ^  ‘ 


(2.8) 


C.  Heat  Energy  Losses 


It  was  shown  in  the  previous  Atlantic  Research  Corporation  report  (i) , 
pp.  4‘-5,  that  the  energy  lost  through  heat  transfer  is  approximately  propor¬ 
tional  to  the  kinetic  energy  of  the  accelerated  mass,  i.  e. , 

[1/2)^-  v2  (2.9) 

where  a  typical  value  of  A'j  is  0.25.  Substituting  this  expression  for  E/f 
into  equation  (2.6)  yields  the  final  form  of  the  energy  balance  equation 


CyC{Ty  -  T)  -Qf-  +  gx  sin  0  . 


(2. 10) 


III.  REDUCED  BALLISnC  SY£^ 

At  this  point,  the  interior  ballistics- of  a  CAD  firing  are  completely 
described  diu'ing  burning  by  the  following  equations. 


Motion 


PA  = 


=  "J  (iLx. 
gc  \  dt^ 


+  g  sin  0 


(3. 1) 


o 


NWL  REPORT  ND.  1752 


State 


NRT 


(3.2) 


Energy  Balance 


Burning  Rate 


Foim  Function 


C  =  +  AjL® 


(3.4) 


(3.5) 


The  syston  is  said  to  have  a  solution  when  all  the  variables  are  deteir- 
mined  as  functions  of  a  single  parameter  such  as  time,  t,  piston  travel,  x, 
or  piston  velocity,  v.  Itofortunately,  the  system  as  it  stands  has  no  solu¬ 
tion  in  closed  fom.  However,  with  the  aid  of  the  following  simplifications, 
the  equations  may  be  reduced  to  a  system  for  which  an  analytic  solution  may 
be  foimd. 


1.  The  nonlinear  burning  rate  eqviation  (3.4)  is  replaced  with  the  lin¬ 
ear  fom 

r  -  OL  Bp  (3. 6) 

where  a  is  chosen  so  as  to  give  the  "best"  linear  approximation  to  equa^ 
tion  (3.4).  One  way  of  achieving  this  end  is  to  choose  a  such  that  the 
integral 

[ap  -  p”)  2  dp 

is  a  minimtm.  This  selection  of  a  ensures  a  "best"  fit  in  the  least  squares 
sense  to  the  burning  rate  over  the  CAD  pressure  range. 
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Carrying  out  the  indicated  integration  of  the  above  equation,  setting 


dcr, 

2^  .=  0,  and  solving  fdr  a  yields 


n-l 


_  3  . 

n  +  2  *inax 


(3. 7) 


where  p  is  the  maxijnum  pressure, 
max 

i^ther  way  of  achieving  a  '^est”  linear  fit  is  to  choose  a  stich  that 
the  integral 


^2  - 


I 


max 


i?|a^  -  I  dp 


da¬ 


is  a  minimum.  Integrating,  setting  =  o  and  solving  for  a  yields 


(3.8) 


Still  another  way  of  obtaining  a  "best”  fit  is  to  choose  a  such  that 

P  P  •  ■ 

/'^max  r  max. 

BP”  dp  -j  •  aBpdp 

0  .  o 

in  other  words,  to  make  the  areas  under  the  two  burning  ratecxirves  equal. 
Carrying  out  the  integrations  and  solving  for  a  yields 

Other  "best"  fits  may  be  derived  in  a  similar  manner.  The  final 
decision  as  to  idiich  definition  of  a  should  be  used  is  rather,  arbitrary. 

Of  the  several  possible  definitions,  equation  (3.9)  seems  tp  lead  to  the 
best  agreement  between  experimental  and  calculated  values  of  maximum  pres¬ 
sure  and  ejection  velocity.  Consequently,  it  will  be  taken  as  the  defini¬ 
tion  of  0. 

Nomographs  depicting  the  relationship  between  a,  p  ,  and  n  expressed 

in  equation  (3.9)  are  given  in  Figures  la  and  lb  in  J^pendix  F."  A  straight 
line  connecting  the  given  values  of  any  two  of  the  parameters  will  pass 
through  the  corresponding  value  of  the  third. 
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2.  6  is  assuned  to  be  0  in  equations  (3.  l)  aod  (3. 3) .  A  suitable  cor>- 

rection  factor  is  developed  in  Section  VI  to  condensate  for  errors  introduced 
by  the  assvndtion. 

3«  The  covolizne  factor,  7)^  is  assumed  to  equal  the  specific  volume  of 

the  propellant.  Thus,  tj  -  ^  =  0  in  equation  (3.2).  Generally,  for  those 

CAD's  operating  under  presstires  below  10,000  psi  the  assumption  leads  to  insig¬ 
nificant  errors.  As  a  matter  of  interest,  however,  a  correction  factor  is 
derived  in  Section  VI  and  may  be  employed  if  desired. 

4.  The  cubic  fom  function,  equation  (3.5)  is  replaced,  for  the  moment, 
by  the  linear  relationship 

C  =  KL.  (3.  ID) 

(Solutions  will  be  obtained  in  Section  X  assuning  a  quadratic  fom  function.) 
This  is  equivalent  to  replacing  the  true  burning  surface  at  any  time  by  the 
average  burning  surface,  S,  during  the  burning  regime.  Thus, 

\  =  pS  -  Cfc/Lfc  .  (3.11) 

For  constant  surface  grains  the  relationship,  equation  (3.  30),  is  of  course 
exact.  ■  Surprisingly,  it  often  leads  to  quite  accurate  results  in  the  case  of 
progressive  and  regressive  grains  as  well. 

It  should  be  noted  that  conditions  2., .3.,  and  4;  need  not  be  included  in 
the  list  of  simplifications;  pn  analytic  solution  may  be  obtained  merely  by 
enploying  condition  I.  The  error,  introduced  by  the  conditions,  however,  is 
more  than  offset  by  the  sinplicity  of  the  resulting  solution.  Any  useful  bal¬ 
listic  system  must  of  necessity  be  a  ccnprcmise'  system  in  which  some  of  the 
accuracy  is  sacrificed  in  order  to  permit  more  rapid  solution  of  problems.  '  '  • 
The  system  developed  here  is  no  exception  to  the  rule.  ’• 

Inposing  the  foregoing  conditions,  and  for  conciseness  writing  ^  ~  • 

Sc 

the  reduced  system  of  equations  to  be  solved  is: 
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During  Burning 


After  Burnt 


PA  .=  V  =  Vv^ 

dt^  a* 

(3. 12) 

(3.13) 

V(r.-n  = 

(3. 14) 

r^§  =  aBP 

(3. 15) 

C  =  kL 

(3. 16) 

(3.17) 

r  12  A*  *  Uil 

(3.18) 

(3.19) 

IV.  DIMEWSICKLESS  SOLUTICTJ:  C  =  \L 


A.  During  Burning 

Conbining  equations  (3.12)  and  (3.15) recalling  that 
dx 

V  .=  2^,  -we  obtain  through  integration;  '  .  :  - 


(4: 1) 


Thus,  the  distance  b\imed  through  the  gj'ain  at  soqy  t.inie  Is  proportional  to 
the*Telocity  of  the  mass  in  motion.'  .Jli  is  here-  assixned •  that  the  mass  begins 
accelerating  iinnediately  upon*  ignition  of  the  propellant,  i.  e. ,  the  value  of. 


^  at  t  =  0  is  zero.  In  view  of  the  fact  that  CAD's  are  essentially  smooth¬ 
bore  guns,  the  assvniiptibn  seems  to  be  a  valid  one,  at  least  for  horizontal 
firings.  Furthermore,  the  a8s\jnf>tion  of  a  non-zero  starting  pressure  overly 
complicates  the  solution;  it  has  therefore  not  been  included  in  the  snalysis. 

Equation  (4.1)  permits  the  reduction  of  the  system  (3.12)  through  (3.16) 
to  the  single  first  order  linear  differential  equation 

(12  +  i/j)  ~  =12  a  BFk  -  6  j  {y  -  l)  A  v  (4;  2) 

up 

where  -y  s  i  +  __  is  the  "ratio  of  specific  heats''^  and  F  =  NRTy  is  the  "in^je- 
tus"‘  of  the  propellant. 

Considerable  sinplification  results  if  we  define  an  adjusted  ratio  of 
specific  heats, 


T  .  1  .  .  1) 


a  dimensionless  distance  variable, 


^  Ui 


and  a  dimensionless  velocity  variable 


rr  —  ^ 

^  "  a'Wk  '' 


for  then  equation  (4;2)  may  be  written 


(1  +  J)  aj  =  1  -  Y  .  (4.6) 

Separating  variables  in  the  equation  and  integrAting,  remembering  that  ini¬ 
tially  X  =  V  =0,  J  is  obtained  as  a  function  of  Y 


( — --i — )  Y- ' 
Vl  .  7.Z.I  y  I 
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The  equation  relates  the  piston  travel  to  its  velocity  at  any  tine  during 
burning.  Figures  2a  and  2b  contain  grsphs  of  I  versus  F  for  various  values 
of  y .  (The  reason  for  incliiding  values  of  y  <  1  will  be  ejqilained  in  Sec¬ 
tion  X. ) 

To  detennine  the  pressure  as  a  function  of  the  piston  velocity,  let  a 
dimensionless  pressure  variable,  F,  be  defined  as 


^  22»\clBF  kj  ^  ‘ 


Substituting  equations  (4.6),  (4; 7),  and  (4.8)  into  equation  (3.12)  then 
yields 


Bof^loying  the  foregoing  equation,  P  versus  7  is  plotted  in  Figure  3  for 
various  values  of  y . 

A  glance  at  Figure  3  will  show  that  P  goes  through  a  mathenaatical  maxi- 
nnm  at  sane  value  V  ■-  Vff.  By  differentiating  P  with  respect  to  F  in  equa¬ 
tion  (4;  9),  setting  the  derivative  equal  to  zero  and  solving  for  F,  the  value 
of  7ff  is  found  to  be 


(4;  10) 


Hence,  the  value  of  P  at  V/f  is 


L 


^JLi 


■  (4.U) 

To  detennine  whether  is  actually  attained  in  a  particular  CAD  firing  we 
must  first  know  the  values  of  two  parameters:  Fj,,  the  'value  of  F  at  grain  . 
burnout,  and  F,|,|,  a  tentative  value  of  F  at  conpletion  of  the  piston  stroke. 


The  valxie  of  F^  is  found  fron  equations  (4:  l)  and  (4.5)  to  be 


b  rvv  * 


(4:12) 
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The  tentative  dimensionless  ejection  velocity  is  obtained  by  a  rearrangement 
of  equation  (4‘.  7)  as 


'ml 


y  -  1 


'(r+ms) 

Hi 


.2LI.I 
2 


(4.13) 


where  is  the  strt^e  length.  (Fj„j  is  a  tentative  value  since  equation  (4: 13) 
is  valid  only  if  burning  takes  place  throughout  piston  travel.) 


The  largest  value  of  P,  call  it  attained  in  a  particular  firing  is 


then 


y  +  1 


p  -  V  / 1  -  7  ^ 

V  y  2  Vj 


(4;  14) 


where  is  the  smallest  of  the  three  values,  and  Vf).  Stated  more 

concisely 


F^  =  min  (Vff,  Ffc) 


(4;  15) 


The  tai5)erature,  T,  may  also  be  detennined  as  a  function  of  the  velocity. 
Substituting  the  values  of  the  dimensionless  variables  into  equation  (3.13) 
with  subsequent  rearrangement  yields 


!•  =  (i  -  k)  . 

B.  After  Burnt 


(4. 16) 


Combining  equations  (3.17),  (3.18),  and  (3.19),  the  differential  equation 
governing  the  ballistics  after  grain  burnout  is  found  to  be 

.  (4. 17) 
(4.4)  and 

(1  +  J)  F  5I  .=  Ffc  -  f2  .  (4:18) 


,  ,  dv  32  FACi,  fl  +  K,  +  K9  \  ,  ,, 

{12  Ax  +  Ui)  V  ~  - -  6A  I - (r  -  1)  • 

Substituting  the  dimensionless  variables  defined  by  equations  (4.3), 
(4:5)  the  governing  equation  sinplifies  to 
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Separating  the  vari^les  and  integrating  from  conditions  at  grain  burnout » 
the  dimensionleaB  velocity  after  burnt  is  obtained  in  the  fom 


V  .=  — 


7-1 


Now,  let 


(4:80) 


i/'X  =  4>  . 


(4.81) 


Then  equation  (4. 19)  may  be  used  to  determine  the  dimensionless  ejection 
velocity, 


=  -J-l  \/%--  *'<■[»  -(i  H--  ’'»)  *] 


(4:88) 


provided  >  Fj,.  In  Figvire  4  are  plotted  values  of  versus  v  for  various 

—  y  —  1 

values  of  y\  Figure  5  contains  graphs  relating  the  quantities  —  ym2> 

The  dimensionless  pressure,  P,  and  the  gas  temperature,  T,  after  pro¬ 
pellant  burnout  are  found  by  algebraic  manipulation  of  equations  (3.17),  (3.18), 
and  (3.19)  to  be 


(4:83) 


r  -  (i  -  |--J  . 


(4:84) 
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Except  for  the  time  ▼ariable,  t,  idiich  capiiot  be  found  in  tenns  of  ele¬ 
mentary  fuDctictns,  *  the  solution  to  the  ballistic  equations  (3.12)  through 
(3.19)  is  DOW  conplete.  The  solution  equations  are  collected  here  for  ease 
of  reference. 

Define  the  dimensionless  quantities 

y  ^  ^ Ar  -  (5.1) 


r  -  12  A 

(5.2) 

^  a.  BF 

(5.3) 

and 

^  12lir  \aiBF  kj  ^  ' 

(5.4) 

The  solution  during  burning  is 

(5.5) 

(5. 6) 

T--Ty  f) 

(5.7) 

/  -  \  7  +  1 

maximum  P  =  (l  -  ^-g--  F^j  7  -  1 

(5.8) 

1  Actually  the  system  (3. 12)  -  (3. 16)  leads  to  infinite  t  in  the  case  P(0)  -  0. 
See  Hirschfelder,  Reference  (b),  page  15,  for  a  discussion  of  the  problan. 
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(5.9) 


(5.10) 


(5.11) 


If  >  F„|,  burning  takes  place  during  the  entire  length  of  the  piston 
stroke  and  no  additional  equations  are  required.  If  Fj,  <  F^j,  buraout  occurs 
before  the  end  of  the  stroke.  The  following  equations  are  applicable  during 
the  period  after  burnt. 


F  = --f-jx/ipT  1-  n)  '  (5.12) 


VI.  .  dOBraiCTIOW  FACTORS  FOR  NOWHORizQMTAL  FIRIMb'aND  OOTOUia  KggBCTS 


A.  Nbnhorizontal  Firing 


Ibe  preceding  analysis  is  based  on  the  assuiDption  that  d  .=  0  in  the 
equations  of  notion  and  energy  balance.  This  is  of  coxirse  not  always  the 
case.  6  may  vary  over  the  range  -!90°  <  6  <  -*90®.  The  magnitude  of  the 
error  introduced  by  the  assunption  is  not  lazge>  being  rarely  more  than  3 
or  4'  per  cent  in  the  most  extreme  cases,  but  a  fairly  reliable  correction 
factor  may  be  derived  and  is  done  so  in  this  section. 


16 


Consider  the  system 


i>A  V  +  g  sin  0^ 
flS  A*  + 

p  I — c — "J  ^  ^ 

CyC{Ty  -  T)  =  I  +  V  gx  Bin  e 


21  -O-BP  (6.4) 

C.=  kL.  (6.5) 

JiB  in  the  case  of  the  friction  force,  Fy,  it  is  assumed  that,  in  the 
large,  the  ballistic  effects  of  the  term,  g  sin  6,  may  be  simulated  by 
replacing  the  term  with  one  that  is  proportional  to  the  pressure,  p.  The 
assim^tion  permits  us  to  write 

^  sin  0  (6.6) 

where  is  a  small  constant  to  be  evaluated  later. 

Substituting  equation  (6.6)  into  equation  (6.1)  and  rearranging  yields 


=  (1  +  A,)  ¥  V  ~  . 

Also,  the  term  gx  sin  6  may  be  written  with  the  aid  of  equation  (6.6) 


gx  sin 


X  X 

in  0  .=  \  g  sin  9  dx  .=  j  ---P-- 

I  .  i  (1  + 


Kj  pAdx 

i  'd  +  YJw 


Substituting  equation  (6.7)  into  equation  (6.8)  and  integrating  we' obtain 

gx  sin  0  •  .  (6. 
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idiich  when  combined  with  equation  (6.3)  produces 

1  +  JT,  +  A", 

- A - i  +  Jf 

Cj,C(Ty  -  T)  =  - -  (1  +  JTj)  I  v2  . 

However,  since  jTj  is  small,  the  foregoing  equation  may  be  very  closely  approxi¬ 
mated  by 

CyCiTy  -  T)  =  (1  +  A'3)  I  v2  .  (6.10) 

A  con5)arison  of  equations  (6.7)  and  (6.IO)  indicates  that  the  net  re.siilt 
of  the  asstxnption  is  to  reduce  the  system  (6.  l)-(6. 5)  to  the  system  (3.12)- 
(3. 16)  with  (1  +  Jrj)^  replacing  V.  Experience  indicates  that  if  ^  ^  sin 

is  less  than  5  or  6  per  cent  of  t  the  assumption  is  a  valid  one. 

max 

Ute  seek  now  the  ballistic  effects  of  a  small  change  in  the  mass  Con¬ 
sider  two  CAD  firii5gs  identical  in  all  respects  except  that  an  effective  mass 
is  accelerated  in  the  first,  and  a  mass  V2  -  (l  +  ^3)li'l  is  accelerated  in 
the  second.  Assuming  that  burning  takes  place  during  the  greater  part  of  the 
stroke,  it  follows  frcxn  equations  (5.2),  (5.5),  and  (5.6)  that  at  any  given 
value  of  the  piston  displacement,  x,  the  quantities  P  and  V  are  the  same  for 
both  firings.  Thus,  it  is  concluded  from  equation  (5.4)  that  at  the  given 
value  of  X 


.Jj__  =  JjL.. 


(6. 11) 


where  subscripts  1  and  2  refer  to  firing  numbers  one  and  two,  respectively. 

2  n  1 

Since  by  definition  a  =  above  may  be  arranged  to  yield 


-  Pi  (1  +  h) 


=  Pi  (1  +  A's) 


(6. 12) 


18 


X  aUi  XYQg 


Tha  equation  must  hold  for  #>j  =  ^  and  •=  Pj^  hence  making  tto  sub¬ 


stitutions  and  solviiK  for  p  in  hav^e 

•  infix  2 


or  since  Jg  <  <  1,  using  teo  tenns  of  the  binomial  e^ansion  gives 


(6. 13) 


(6.14) 


The  effect  then  of  multiplying  the  mass  V  by  a  factor  (1  +  fg)  is  to  multiply 


the  maTimin  pressure 


eqiproximately  by  a  factor  f-A) 


To  detennine  velocity  effects  we  note  that  equation  (5.3)  permits  us  to 


write 


'Ll  ^ 

“l  “2 


(6. 15) 


at  any  given  value  of  x.  Using  the  definition  of  o,  the  above  then  yields 


1  '^maxi/ 


I 

3  -  2n 


V 1  •  (6. 16) 


Ihe  relation  must  hold  at  x  -  x^,  hence  again  using  two  terms  of  the  binomial 
expansion  we  have 


''m2  « 


^1  -  Sj-i-g-  v„j  .  (6.17) 


(As  a  rather  Interest  i^  by-product  of  the -analysis  it  is  noted  that; 

1.  For  a  biuming  rate  exponent  of  one,  the  ejection  velocity  is  inde¬ 
pendent  of  changes  in  the  accelerated  mass; 

2.  For  exponents  greater  than  one,  the  velocity  increases  with  increas¬ 
ing  Vi 
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3.  Exponents  much  largM*  then  one  should  be  aroided  because  of  the 
extreme  pressure  sensitivity  as  n  1.5. ) 

The  vslue  of  JTj  to  be  used  to  conpensate  for  nonhorizontal  firing  is 
obtained  from  equation  (6.9)  by  setting  v  =  *  ~  *m' 

f.il'j.JL*?./ »  ,  (6.ie) 

"m  »i.i 

B.  Covolume  Effects 

In  the  lazge,  the  ballistics  effects  of  the  covolume  term  may  be  simu¬ 
lated  in  the  following  manner.  Consider  the  equation  of  state  (3.2)  in  the 
form 

P  -  M-n  -  ^)  =  12  W  .  (6.19) 

p 

If  now  the  tern  P(t)  -  — )  is  replaced  by  an  average  value  (tj  -  -^) , 

Instead  of  assixnlng  v  ~  ^  ~  0,  equations  (3. 12)  -  (3. 16)  reduce  to  the  single 
differential  eqiiation 

(12  iU  +  tfj)  ^  .=  12  a  BF '  X  -  6  (r  -  l)  ,  (6.20) 

where 

*  W  7>] 

Equation  (6.20)  is  identical  to  equation  (4.2)  with  the  exception  that 

p 

F'  replaces  F.  Furthermore,  the  term  (■>?  -  ■^).  is  generally  less  than 

0.02  for  most  CAD's.  Thus,  the  error  introdticed  by  neglecting  the  covolume 
term  may  be  approximately  coopensarted  for  by  a  sm^l  increase  in  the  impetus, 
F,  of  the  propellant. 

To  determine  the  ballistic  effects  of  a  small  increase  in  F,  again  con¬ 
sider  two  CAD  firiiigs  having  identical  ballistic  parameters  except  that  the 
inpetus  of  the  first  is  Fj  and  that  of  the  second  is  F2  -  +  K^)F^,  where 

K^-  <  <  1.  As  in  the  preceding  analysis,  it  is  concluded  that  -  at  any  given 
value  of  the  di^laoeroent,  x, 
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F|  *  f  2  *  ttj* 


(6.22) 


where  again  aubscripta  1  and  2  refer  to  firing  nuaijers  one  aitid  two,  respec- 
tiTely.  Bifiloying  the  definition  a  =  relationship  between 

the  two  maTiwin  pressures  is  found  to  be 

_  2 

.  9  . 

P 


max  2 


=  (1  +  ^4)*“*"  P 


max  1 


(6.23) 


or,  approximately. 


^max  2  ^  Z-Zn/^Toaxi 


(6.24) 


To  detennine  the  effect  on  the  ejection  velocity  it  follows  from  equa¬ 
tion  (5.3)  that 


=  .ymx. 


a,F 


l*’! 


«2^2 


(6.25) 


so  that  using  equations  (6.22)  and  (6.23) 


^"2  -  (1  +  ^ 
V  ^max  1 


)»-2n^ 


ml 


Hence, 


^m2  *  ^  3~-^27i)  ‘  • 


(6.26) 


(6.27) 


The  value  of  JT^  to  be  used  to  ccnpensate  for  covolume  effects  is  detei^ 
mined  from  equation  (6.21).  '  ^  . 


r  ■  =  _max_l  (t,  _  i)  «  -Wl  (ji  _  . 

*  24F  .  o'  ^  .  '24^  o' 


(6.28) 


C.  Combined  Correction  Factor 


The  corrections  for  nonhorizontal  firing  and  neglected  covolume  may  be 
combined  into  a.  single  correction  factor.  The  following  procedure  is  suggested. 
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First  assume  6-0  and  tj  =  — .  Calculate  the  maxlimiti.  pressure, 
and  ejection  velocity,  using  the  method  described  in  the  foUoHiog  section. 
The  corrected  maxinun  pressure,  and  the  corrected  ejection  velocity,  v*, 

are  determined  from 

=  (l  ♦  (l  *5V-E,-  <«•=«> 

and 


where  .  . 


2  sin  9 


and 


VII.  (aiAHgCAL  EETEFMINATICK  OF  VJiXBim  PRESSURE  AND  EJECTION  VELOaTY;  C  =  \L 


The  fact  that  the  parameter,  a,  is  itself  a  function  of  the  maximum  pres¬ 
sure  precludes  an  explicit  determination  of  the  maximum  pressure  and  ejection 
velocity.  However,  a  rapidly  converging  trial-and-error  procedure  may  be 
devised  utilizing  the  equations  summarized  in  Section  V.  The  procedure  requires 
more  than  one  iteration  only  if  the  maximum  pressure  occurs  at  ’Taumt". 

We  begin  by  defining  the  parameters: 


a  -  M  r 

12¥  \BFJ^J 

(7.1) 

„  -  A 
^2  bF\ 

(7.2) 

a  - 

^3  F\M 

(7.3) 
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1# 

Cl  “  «7i/a* 

(7.4) 

atwl 

0,  »  qj/a  . 

(7.5) 

The  dimensionless  vsriables  defined  by  equations  (5.3),  (5.4),  and 
(5.11)  are  then  written  simply 


P  =  QiP  =  <7i^/a* 

(7.6) 

V  ■-  CjV  =  q^v/a 

(7.7) 

and 

(7.8) 

Now,  p  may  be  eliminated  from  the  two  equations 

max 

'’^  =  ’1  W* 

and 

_ 2  .  n-i 

(X  — T-t  P 

n  +  1  '^max 

to  yield  a  as  a  function  of  the  burning  rate  exponent,  n. 

Thus 

and  the  ratio 

(7.9) 

The  relationship  between  a,  n,  and  expressed  in  the  foregoing  equation 

la  shown  nomographically  in  Figures  6a  and  6b.  (Again,  a  straight  line  con¬ 
necting  the  values  of  any  two  of  the  parameters  passes  through  the  corre¬ 
sponding  value  of  the  third.) 

The  trial-and-error  procedure  may  now  be  outlined. 
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step  1 


a.  Usipg  the  given  data  for  a  particular  firing,  calculate  the  values 
of  the  parameters 

X.  =  Cj/Lj 

y  =  (y  .  1)  ,1 

Ui  =  Vc  -  Ck/  P 
V  =  (1  +  K^)ur/g^ 

,  =  Hi. 

12lir  \BFkl 

n  =  * 

^2  Wk 


a.  Assume  a  valvie  of  P,,.  (initially,  this  is  done  by  selecting  the 

^  • 

value  of  P  fixm  Figure  3  corresponding  to  the  smaller  of  the*  two  values  yml 
and  'Fj^as  determined  in  Step  1. )  '  .  ‘ 
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b.  Calculate 

c.  Detemizie  a  flxxn  Figure  6  corresponding  to  the  bumii^-rate  exponent, 
n,  and  the  ratio,  P.^/q^. 

d.  Calculate =<73  /0*. 

e.  Determine  =  min  (Fj^  Fj,  F«l).  Select  the  value  of  =  /(F^  ^ 
from  Figure  3. 

f.  :if  the  assuned  value  of  P^  equals  the  calculated  value,  go  on  to 
Step  -3;  if  not  repeat  Step  2  using  the  calculated  value  of  P^ 

Step  -3 

a.  Calculate 

Ql  =  <7i/a* 


and 

(?j  =  <7j/0  • 

b.  Calculate  the  maximum  pressxire  using  ~  Pf^/ Q\' 

Step  4 

Conpare  the  values  of  Vf^  and  ’Rel¬ 
ease  A;  Vf,  Calculate  the  ejection  velocity  using  •=  WC2- 

Case  B:  Fj,  <  F„j. 

a.  Detemine  Jj,  .=  /(F5,  y)  from  Figvire  2. 

1  +  Jfc 

b.  Calculate  v  =  • 

c.  Determine  ^  =  /(v,  7)  from  Figure  4'. 

d.  Determine  F„,2  =  ^  Figure  5. 
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e.  Calculate  F,,  =  . 

f.  Calculate  the  ejection  velocity  using  .=  Fj,2  /(?2  • 

Step  5 

The  correction  factors  derived  in  Section  VI  may  be  applied  if  desired. 
Calculate 

.  r  =2_Fx„'8in_0 

*3  “  '•“~2 

and 

-  The  corrected  values  of  maxiimm  pressure  and  ejection  velocity  are  then 
detemined  from 

*  rrrn (l  ♦  5“Vn  ^«) 

and 

■  3””2>j  ^3)  3"=^"2n  * 

San5)le  Problem  1. 

Given- the  following  data 
i4.=  4;92in2 

5=  0.02101  in/sec 

Cfc  =  0.161  lb„^ 


includes  the  main  charge  weight  pl\is  the  "equivalent  .igniter  weight."  The 
latter,  denoted  by  Cj,  is  defined  as  Cj  .=  [F^/ F)l,  where  F,-  is  the  impetus  of 
the  igniter  charge  and  J  .is  the  igniter  charge  weight. 
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F  =  205,500  ftaby-Zlb^ 

1 1  =  0.  ID 
Jj  =0.25 
Li,  =0.075  in 
n  .=  0. 47 
Oc  =  159. 1  in? 
w  =  360  Ib^ 

=  2.725  ft 

0  =  70. 5  degrees 
y  =  L235 
p  =0.06  lb„/in? 

77  =-30  in3  /lb„ 

determine  the  maxiimxn  pressure  and  ejection  velocity. 

The  given  data  are  listed  in  Table  1.  The  required  parameters  are  then 
calculated  and  the  values  Fj^  •=  0. 7764' and  F^j  =  0.6800  are  obtained.  The 
initial  assxmed  value  of  is  taken  to  be  =  0.295.  The  assumed  and 

caloulated  values  of  are  found*  to  be  equal  after  one  iteration  and  the 
maxisnxn  pressure,  is  calculated  to  be  1496  psi.  Since  Fj  >  F„j,  the 

ejection  velocity  is  V^i/Q^  ~  49. 8’ ft/sec.  The  correction  factors  of  Sec- 
tion  VI  are  then  applied  leading  to  =  3£6Q  psi  and  .=  49.0  ft/sec. 

The  values  conpare  favorably  to  those  obtained  through  numerical  integra¬ 
tion  of  Equations  (3. 12) -*(3.16) -namely,  1607  psi  and  51.9  ft/ sec. 

i\n  estimate  of  the  errors  inherent  in  the  graphical  method  can  be  made. 

Assvming  and  a  are  determinable  from  the  graphs  to  within  i  2  per  . 

cent  of  their  correct  values,  the  maximum  errors  in  p  and  i*  are  about 

.  ’  max  ® 
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i  6  per  cent  and  t  4*  per  cent  req)ectiTely.  This  obsenration  followa  traa 

P  Y 

the  act  that  p  =  -^  o?  and  v_  .=  ~i  a. 

''max  ■  7j 

San^le  Problem  2* 

Given  the  folloiring  data 

A  =  O.^asd-ii^ 

B  =0.01741  (ii^/lby)"  In/aec  .  ‘ 

Cfc  =  0.009304  lb*  ■  •  ■ 

F  =  362,300  ft -Iby/lb* 

=0.1  ■  .  • 

=  0.25 

Lfc  =  0.055  in  ■  ,  •  . 

n  =  0.51 
Uf.  =  6  in’ 
w  =  2100  lb„ 

=  1  ft 

y  =  1.215 
0=0" 

p  =  0.06146  lb* /in’ 

77  =  30  in’ /lb* 

detennine  the  maxiimmi  pressure  and  ejection  velocity. 
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The  given  data  are  listed  in  Table  2.  ;in  this  case,  the  maximum  pressure 

occurs  at  burnt  so  that  five  iterations  are  required  to  obtain  convergence  of 

the  P,,  values.  The  results,  P  =  6367  psi  and  V-  .=  8.7  ft/sec  again  com- 
M  max  >“  / 

pare  favorably  with  the  values  6619  psi  and  8.9  ft/sec  obtained  from  numerical 

integration. 


VIII.  DETERMINATION  OF  CHARGEE  WEIGHT  AND 
Wm  SIZE  TO  MEET  GIVIN  PERFORMANCE  SFECIEI CATIONS 


A.  Given  p  and  v- 
max 

If  it  is  specified  that: 

a.  The  propellant  shall  bum  throughout  the  piston  stroke,  i.  e. , 

and 

b.  p  shall  be  a  mathematical  maximum,  i.  e. ,  V„y  >  Vu  =  -4  » 

max  ^  y 

the  charge  weight  and  web  size  required  to  produce  a  given  maximum  pressure 
and  ejection  velocity  may  usually  be  calculated  explicitly.  Conditions  a. 
and  b.  are  not,  ir.  general,  serious  restrictions.  They  are  usually  called 
for  in  roost  CAD  applicatipns. 

Let  U(.  be  the  volime  of  the  entity  chamber.  By  definition  of  initial 
free  volume  we  then. have,  using  equation  (5.2) 

i/i  =  (8.1) 

and  solving  for 

Cj  -  y.  .  (8.  2) 

Condition  a.  above  permits  equation  (5.3)  to  be  rearranged  to  give 

^  ^  "4  " 
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'ml 


(8.3) 
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or 


h 


.(oiBF\  Vnj 


Cb 


(8.4) 


1^,  is  related  to  expression 


Xm- 


-A. 

vy- 1 


1  -  ^Ll-l  Y 


-  1 


(8.5) 


hence,  if  can  be  determined  frcm  the  given  conditions,  equations  (8.2)  and 

(6.4)  determine  the  required  charge  weight  and  web  size.  We  proceed,  then,  to 
determine  as  a  function  of  the  maximun  pressure  and  ejection  velocity. 

Combining  equations  (7.1),  (7.2),  (7.4),  and  (7.5)  we  obtain  the  rela¬ 


tion 


i2"F  ^2 


(8.6) 


Uultiplying  nixnerator  and  denominator  of  the  right-hand  member  by  vj,  and 
using  equation  (7.7),  this  may  be  written 

0  . 

*  12  *  (8.7) 

By  condition  b.  above,  equations  (4.11)  and  (7.6)  may  be  combined  yield- 

7  U 


ing 


^  >nax  y\Zy) 


(8.8) 


Combining  equations  (8.7)  and  (8.8)  and  substituting  for  its  equivalent, 
12  ^^1 1^1  we  obtain 

A 


2  +  1 


^roax  _  _1  /f  ~  ^ 

y  \  ”2’^' ) 


(8.9) 


which,  when  solved  for  X  yields 


Xii  r  1 

r  =  ^  y  -  ^  \  ^max 

"  M2  +  i)  [(i/2)Wv„^J 


"ml 


(8.10) 
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Define  tl»  two  paraoieterSf 


and 


^ 


(1/2)  Vp* 


(8.11) 


(8.12) 


^Piezooetric . efficiency  is  defined  by  comer  (d)  as  the  ratio  of  the  mean  pres¬ 
sure  to  the  mwTlninn  pressure.  'Hie  mean  pressure  . in  this  definition  is  the 
pressure  idiich,  when  applied  to  the  accelerated  mass  over  the  total  strcdce  in 
the  device,  will  produce  the  observed  velocity.  Hie  maximum  pressure  is  the 
observed  maxiiami  pressure. 

Equation  (8.12)  may  be  derived  as  follows: 


Hie  kinetic  energy,  of  the  mass  at  the  . end  of  stroke  is 

Bit  =  1/2  Vv* 

and,  by  the  above  definition  of  mean  pressure,  p, 

Ek-  • 

Equating  the  two  expressions  for  and  solving  for  f, 


-  W-v- 

P  = 


Dividing  by  p 


max 


IT  =  -t.  =  ...lyji,. 

^max 


max 


H»  ratio  p/p  gives  an  indication  of  the  flatness  of  the  pressure- 

travel  curve.  In  a  device  with.no  retarding  forces  (a  completely  friction¬ 
less  device  accelerating  a  mass  horizontally),  the  ratio  would  give  a  true 
indication  of  the  flatness  of  the  pressure- travel  curve.  'Pb  retarding 
fomes  increase  in  the  device,  the  ratio  indicates  less  certainly  the  shape 
of  the  pressufe-travel  curve.  In  the  case,  for  exanple,  where  the  retard¬ 
ing  forces  lequaled  the  force  accelerating  the  mass,  a  piezometric  efficiency 
of  only  50%  would  be  obtained  from  a  pressure-displacement  curve  having  a 
joenstant  pressure  tbrou^mt  the  stroke. 
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Equation  (8.10)  then  reduces  to 


r  =  JC'7  * 

g  'ml 


(8.13) 


But  is  also  related  to  by  equation  (6.5).  Ifence,  for  any  given  value 
of  E,  7^^  can  be  determined  frcm  the  liq>licit  relationship 


Ly  »  =  / 

E  r 


1  -  7 

2 


y  - 1 


(8.14) 


Xjn,  Cff,  and  Lf,  may  then  be  obtained  from  equations  (8.5),  (8.2),  and  (8.4), 
respectively. 

There  aire  two  real  roots  of  equation  (8.14)  corresponding  to  each  value 
of  E  and  y;  however,  the  larger  of  the  two  roots  generally  leads  to  prohibi¬ 
tively  large  values  of  C{,  and  Lj,  and  may  be  disregarded,  i^pendix  E  treats 

the  cases  in  which  the  larger  root  must  be  tised.  In  Figure  7a,  we  have 
plotted  the  smaller  root,  which  is  denoted  by  Xnl  versus  B  for  several  values 

of  f.  The  curves  labeled  E  and  E  .  are  such  that: 

max  min 

if  E  >  E  the  roots  of  equation  (8.14)  are  imaginary; 

inouK 


if  E  <  E  .  , 
min’ 

b.  is  violated. 


the  smaller  root,  is  less  than  l/y  =  7ff  and  condition 


The  problem  of  determining  Cj,  and  Lf,  to  produce  a  given  value  of 
E  =  / v^)  is  said  to  have  a  ’’feasible  solution"  if  satisfies  con¬ 
ditions  a.  and  b. ,  above,  i.  e. , 

F6  >  >  1/f  . 

The  procedure  for  determining  the" charge  weight  and  web  size  may  now 
be  outlined. 

Step  1 

9^  * 

a.  Given  the  maximum  pressure  and  the  ejection  velocity 

together  with  the  other  ballistic  parameters  for  a  particular  firing, 
calculate  the  values: 
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Step  3 

Calc\ilate  the  values 

7  -  o^BV  ■ 

min  ''m 

Cf,  -  piUf.  -  Ui),  and 
Lb  =  Cfc/X  . 

The  foregoing  values  give  the  required  charge  weight  and  web  size  pro¬ 
vided  Lj,  ■>  Lfj  miHw  •  If  Lfj  <  Lf,  njin  the  problem  has  no  feasible  solution. 

Saople  Problem  3. 

Other  data  being  the  seme  as  in  Semple  Problem  2,  what  web  and  charge 
weight  are  required  to  yield  =  6000  psi ,  and  .=  10. 5  ft/sec?  The 

solution  is  shown  in  Table  3.  It  is  concluded  that  a  constant-s\irface 
charge  weighing  0.0&4  Ib^  with  a  web  of  0.73  in  will  meet  the  specifica¬ 
tions.  Notice,  however,  that  the  design  leaves  quite  a  bit  of  unbumed  pro¬ 
pellant  at  ejection.  The  distance  burned  throxigh  the  grain  at  the  end  of 
the  stroke  is  only  0.30  in,  i.  e. ,  about  41  per  cent  of  the  web.  A  dis- 
cussicm  of  this  problem  is  given  in  Section  IX,  and  i^pendix  E. 

B.  Given 

Frequently,  the  ejection  velocity  is  specified  without  attaching  much 
iiqportance  to  the  maxlimm  pressure  provided,  of  course,  it  stays  within 
certain  safety  limits.  In  such  cases  several  values  of  E  may  be  selected 
in  the  peimissable  range  and  the  corresponding  maximun  pressiire  determined 
from 

^max  ~E  Ax„ 

Values  of  Cj,  and  Lf,  may  then  be  calculated  in  the  usual  way.  This  slight 
revision  of  the  method  has  two  advantages: 

1.  For  any  selected  value  of  E  it  can  be  determined  immediately 
whether  or  not  a  feasible  solution  exists,  and 
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Table  -3  ngmamMTTnw  OF  CISUK2;  WEIOir  AND 

WEB  SIZE  TO  MEET  GIVM  PEBFOEaMNCE  SFECIFrCATIOWS 


SaD;>le  Problem  3. 

Given  Values 

Step  1 

A  =  0.7854  in2 

B  =0.01741  (in*  /lby)"  in/sec 

F  =  3.623*  ft.lby/lb„ 

Ky  -  0. 10 

=  0.25 

«  =0.51 

tfc  =  6  in* 

=  10.50  ft/sec 

u>  =  2100  lb„ 

=  1  ft 

^3=2  gx^  sin  =  0 

P*  (77-^) 

=  0.0095 

P  =  P  *  /(l  +  (1  +  --*)  =5943 

'^max  ^max/  3-2n/  V  3-2«/ 

-n  =  /(l  -  (1  ^3=80=  “-^5 

F  =  (1  +  ^’j)  =  71.797 

fl+K,  +  JT,  \ 

f  =  1  1  (y_i)  +  1  =  1.2639 

£  =  (1/2)  M  v^lP^Ax^  =  0.0396 

y  =  1.215 

Step  2 

6=0  degrees 

a  =  /  ^^max’  ^^Eure  1  =  1.82'2 

sin  6  =  0 

V„Y  =  f  {E,  y)\  Figure  7a  =  0.97 

p  =  0.06146  lb„/in* 

Tj  =  30  in*  /  lb„ 

~  ^  ^%tV  J  Figure  2  =  1.84 

Ui  =12  kx^l  -  5.122 

(Note:  0^  denotes  a  x  10^) 

^  =  0.0737 

Step  3 

^^n,in  =0.  3087  to 

Ci,  =  p  [n^  -  Ui)  =  0.05396  lb„ 

Ih  =  Ch/K  =  0.732  in 
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2.  A  wide  range  of  pemisaable  values  of  thew^  size  and  charge  weight 
may  be  determined  which  yield  very  nearly  the  same  and  v„ . 

Using  the  inequalities 

Ffr  >  r*,  > 


it  can  be  shown  that  for  any  selected  value  of  a  feasible'solution  exists 
if  and  only  if  satisfies  the  relationship 


(l/2)»r_v„2  F^J 

F  p  Uc  -  "2 


(lAm) 


where 


(Note  that  Figures  7a  and  2  may  be  used  to  evaluate  the  right-hand  member 
of  the  inequality;  also  if  satisfies  the  inequality  for  a  particular 
value  of  E  it  will  satisfy  it  for  all  larger  values  of  E  since  the  right- 
hand  member  is  a  monotonically  increasing  fimction  of  E.) 


The  alternative  method  may  be  summarized  as  follows. 


Step  1 

Given  the  ejection  velocity,  v*,  and  other  ballistic  data,  calculate 

.  _  2  sin  9 

"3 - vr - 


=  V- 
In  Ht 


3"-"2n  ^3/ 


=  (1  +  Ki)Wec 


^The  covolune  connection  factor  unduly  complicates  the  method  and  so  is 
not  included. 
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_  (i/2)_»rv„2 

*1  ■  ~rpur~ 


step  2 


a.  Select  a  value  of  £  in  the  range  <  E  <  from  Figure  7a 


b.  Determine 

=  f(E,  y)  from  Figure  7a  and 
~  Figure  g. 

c.  Calculate 


^2 


(1/4) 


If  £j  <  fj  6°  Step  3.  If  not,  select  a  larger  value  of  E  and 

repeat  Step  2. 


Step  5 

a.  Calculate  t 

'^max 


(1/2)  m  v^_ 


b.  Determine  a  .=  fip  ,  n)  from  Figure  1  , 

max 

c.  Calcvilate  the  values 

f/i  =  12  AxJX„ 

^  ~  oTBF 
C{,  =  p{Uc  -  Ui) 

and 


Lb  =  C<,/\  . 
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d.  Rqpeat  Steps  2  and  -3  for  sereral  different  ralues  of  E.  Plot 
versiis  on  one  or  two  cycle  log  paper.  The  resulting  curve  gives  the 
required  charge  weight  as  a  function  of  the  web  size  to  yield  the  given  ejec¬ 
tion  velocity. 

Ttie  results  of  a  typical  calculation  are  plotted  in  Figure  8. 

Notice  the  insensitivity  of  to  large  changes  in  and  L/,.  In  the 

pressure  range  1600  <  p  *  <  1650  (a  variation  of  about  3  per  cent)  the  web 

B)8Z 

and  charge  weight  vary  by  350  per  cent!  This  is  a  phenomenon  comnon  to  most 
CAD's  where  low  loading  densities,  i.  e.,  C^f  «  1,  are  the  rule. 


C.  Given  p 
_ max 

A  method  may  also  be  devised  to  handle  the  case  of  a  given  maximum  pres- 
siire,  with  an  unspecified  ejection  velocity.  The  method  is  only  a  slight 
revision  of  the  one  given  in  part  B  of  this  section. 

Step  1 

Given  the  maximun  pressure,  p  *  ,  and  other  ballistic  data,  calculate 

inax 

^  34  f 


F  =  (1  +  Ki)u,/g^ 


and 


y  = 


(7-1)  +  1  . 


Step  2 

a.  Select  a  value  of  £  in  the  range  E  .  <  E  <  E  from  Figure  7a. 

min  -  max 

b.  Detemine 

=  /  (£,  f)  from  Figure  7a,  and 
^  y)  from  Figure  2. 
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d.  If  <  E2  go  Step  3.  If  not,  select  a  larger  value  of  E  and 

repeat  Step  2. 

Step  3 

a.  Calculate  ~ 

m 

b.  Determine  a  .=  n)  frcoi  Figure  1. 

c.  Calculate  the  values 


>r—[i7^r 


■-  12  Ax^/ X„ 


Ct  =  p{Uc  -  Ui) 


and 


h  =  Cfc/X  . 

d.  Repeat  Steps  2  and  3  for  severs  1  values  of  E-  A  plot  of  Cj,  versus 
on  one  or  two  cycle  log  paper  gives  the  required  charge  weight  versus  web 
size  to  yield  the  given  maximum  pressure. 

IX.  BALLISnC  EFFECTS  OF  SHALL  PARAMETER  CHANGES 

A  useful  outgrowth  of  the  equations  simriarized  in  Section  V  is  the  fact 
that  if  one  or  more  of  the  loading  conditions  of  a  particular  CAD  is  changed, 
the  corresponding  change, in  the  performance  may  be  calculated  and  vice  versa. 
This  fact  is  particrilarly  useful  in  modifying  grain  designs  that  do  not  quite 
meet  performance  spec if icatioos.  For  exan^le,  svqppose  the  specifications  for 
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a  CAD  call  for  a  mat  limn  pressure  of  6000  psi.  Using  the  method  of  Sec- 
ticQ  VIII  a  charge  weight,  C^|,  and  a  web,  £51,  are  detemined  idiich  should 
meet  the  specifications.  A  charge  is  made  up  having  the  indicated  dimensions 
wnH  is  fired  in  the  device.  Ihe  measured  maximum  pressure  is  5000  psi.  Uhat 
changes  in  the  grain  design  should  be  made  to  meet  the  specifications? 

Assuning  that  biixntng  tdces  place  throughout  the  greater  part  of  the 
stroke  length,  and  that  a  mathematical  pressure  maximnn  is  attained,  equa¬ 
tions  (5.8),  (5.9)  and  .  (3.9)  permit  us  to  write 

where  snd  are  the  new  values  of  the  charge  weight  and  web  size 
required  to  attain  6000  psi.  More  generally,  if  it  is  desired  to  change  the 
maximum  pressure,  p  . ,  to  a  new  value  p  .=  (l  +  e)  ^  .  where 

I  €  I  «  1  we  may  write 


(1  +  -  2n 


There  are  then  three  possible  modifications  of  the  charge  design. 
(1)  Keep  the  same  surface  area  and  increase  the  charge  weight.  In  this 


^  hi  h2 


so  that  substitution  in  equation  (9.S)  yields 

-  ^61  =  iPh  -  hz)  (1  +  €)' 


Defining  A,  the  loading  density,  by 


Equation  (9.4)  may  be  solved  for  C{,2  to  yield 


(9.5) 


^62  •“  ^61 


[i  -  (i  - ") 


(1  +  £)-(»  - 
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or  approxlinately 


Cbi  •  +  (»-2n)  ej  Cfc,  .  (9.6) 

nram  equaticn  (9.-3)  «e  also  obtain 

Equations  (9.6)  and  (9.7)  then  determine  the  near  values  of  the  charge  weight 
and  web  size.  (Notice  that  this  modification  is  practical  only  if  A'»  0; 
otherwise  the  equations  may  lead  to  very  large  values  of  ^bt’^ 

(Z)  Keep  the  same  web  and  increase  the  charge  weight. 

In  this  case  we  may  substitute  into  equation  (9.2) 

^bi  •"  ^bi  (9*8) 

to  obtain 

“r^  -  (1  +  6)»  -  *« .  (9. 9) 

.  ^b2 

Solving  for  Cf^  then  gives  epproximately 

(3)  Keep  the  same  charge  weight  and  decrease  the  web. 

In  this  case 

Cb2  =  Cfci  (9.11) 

and  equation  (9. 2)  reduces  approximately  to  , 

Ltt»  •  (9-12) 

In  a  similar  manner,  the  variations  In  maximum  pressure  and  ejection 
velocity  due  to  changes  in  the  other  ballistic  parameters  may  be  determined. 
The  following  table  simnarizes  the  results  of  such  an  analysis. 
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Nm 

!  - 

Table  4-  BALLISnC  YARIATKK3  DOE  TO  SMALL  EAMMEIER  CHANCES 


1  Per  Cent 

Per  Cent  Variation 

Per  Cent  Variation 

Increue 

in 

in 

in; 

Uaxioun  Pressure 

Ejection  Velocity 

A 

f  2  \ 

V3  -  Bn) 

3-2"  1+1 

B 

2 

1 

1 

3  -  Bn 

/  1  \  ^  A  \ 

1  +  /  A  \/  S  l-n\ 

3^&J  U-A/Vi+^  3-^/ 

\3  -  2nH^  1  -  aJ 

e 

2 

1 

3  -  2n 

3  -  2n 

f  , 

/  2  \ 

/  1  \ 

\3  -  2n; 

-  (vh) 

l3  -  2"] 

/  1  \  /  1  -r.  AS  \ 

\3"-"^  ■  r"5; 

1 

/I  -  n  \ 

3  -  2n 

\3  -  2n; 

S 

0 

f+T 

f-  -  •  ■ 

Note; 

O 

II 

and  S  .=  -5" — 

12  Ax„ 

It  should  be  eoq)hasized  once  again  that  the  detemination  of  the  varia¬ 
tion  effects  is  based  on  the  assunptions:  C  =  X.L ,  and  -  *^^6  •  How¬ 

ever,  for  small  changes  in  the  parsmeters  the  indicated  variations  in  maxiimm 
pressure  and  ejection  velocity  should  be  of  the  right  order  of  magnitude  for 
most  CAD  firings. 

Another  useful  relationship  may  be  derived  frcm  equation  (9.2).  If  the 
ratio,  yS,  defined  by 
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is  held  constant,  the  niaxijiium  presstire  is  imchanged  provided  ^  F^. 


Rir- 


thennore  for  low  loading  densities. 


i  e 

pU, 


«  1,  the  ejection  velocity  is 


essentially  unchanged  for  constant  ’/S.  (This  follows  from  equations  (5.2), 
(5.5),  snd  (5.3).)  Tliis  relationship  is  useful  when  it  is  desirable  to  attain 
a  given  maximum  pressure  (or  ejection  velocity)  and  to  have  "burnt”  occur  at 
the  moment  of  coa|)letion  of  the  stroke.  To  achieve  that  end  the  following 
procedure  may  be  used. 


A  value  of  v_  (or  i>  )  is  assumed  that  leads  to  a  value  of  £  in  the 
^  max 

range  E  ,  <  E  <  E  .  U4  and  k  are  calculated  using  one  of  the  methods  in 

min  “  “  max  » 

Section  VIII.  Hie  values 


^£-Cb/P^Ut 

^  W/Lb)^ 


(9. 13) 


htiH 

4min  =("r)’^m  0.14) 

are  calculated,  is  then  the  required  web  size.  Solving  for  Cj,  in 

equation  (9.13)  then  yields  the  required  charge  weight. 


Cb 


(9. 15) 


X.  DIMENSICaJLESS  SOLUTION;  C  =  \,i  + 

Similar  solutions  to  those  described  in  equations  (5.5)  through  (5.14) 
may  be  obtained  if  the  linear  form  function  C  ~  XL  is  replaced  by  a  quad¬ 
ratic  form  function  of  the  type 

C  =  \iL+kzL^  .  (10.1) 

Before  elaborating  further,  let  us  evaluate  the  coefficients  Xj  and  Xj 
in  the  equation.  As  before,  let  X  =  We  then  obtain 

X  =  Xj  +  Xj  Lj, 


45 


IIWL  HEPORT  ND»  1752 


or 


k- X., 

"h  ~  * 


(10.’2) 


Furthermore,  since 


dC 

az 


=  PS 


(laa) 


where  S  is  the  burning  siirface  at  any  instant,  we  also  have 


d*C_  .d5  _ 


dt 


(10.4) 


dS 


Replacing  with  an  average  value. 


(St -  So) 


,  where  S^,  and  St  are  the 


values  of  5  at  I  .=  0  and  L  -  L5,  respectively,  and  substituting  into  equa¬ 
tion  (10.2)  yields 


k^  =k~  piSt  -  So)  /2 


(10.  5) 


Equations  (10.2)  and  (10.5)  are  then  the  values  of  \j  and  Xj  to  be  used  in 
equatlcm  (10.1).  (As  a  matter  of  interest,  it  can  be  shown  using  the  above 
equations,  that  if  the  true  form  function  is  a  cubic  of  the  type 


C  =  A,i  +  AjA*  +  A,i* 


the  values  of  Xj  and  Xj  are  given  by 

X|  -  A|  -  (1/2)  Aj  Lt 


and 

Xj  —  A2  (3/2)  Aj  . ) 

Values  of  Xj  for  some  of  the  more  comnon  grain  shapes  are  given  in  Appendix  A. 
Equation  (lO.'2)  may  then  be  \ised  to  determine  Xj  . 

After  b\imt,  the  ballistics  are  independent  of  the  form  f\jnction,  hence 
we  need  only  concern  ourselves  with  the  period  before  b\mit.  During  the  buuvi- 
Ing  regime  the  system  to  be  solved  is 

PA  =  Vvil  (10.6) 
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(10.7) 

-  n  ■ 

i  »’ 

(10.8) 

(10.9) 

and 

C  =  X,  I  +  XjL*  . 

(10.10) 

Oonbiniog  equatioos  (10.6)  and  (lO.  9)  we  again  obtain  the  relation 

L  =  V  .  (10. 11) 


Ihe  system  may  then  be  reduced  to  the  single  differential  equation 
(12  Xx  +  ffj)  ^  -  12  aBf  +  12  (aB)2  v 

-  6  (y  -  1)  Av.  (10.12) 

Defining  y  and  I  as  before, 

r  =  1  +  (y  -  1)  (10.13) 

and 

I  =  ---  X  (10. 14) 

but  with  a  new  velocity  variable,  F',  defined  as 

■  aix-  «  ■  (10- 15) 

equation  (lO.  12)  simplifies  to 

(1  +  I)  =  f  •  (10-16) 
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Udir,  eoployiog  equations  (10.11)  and  (10.15)  t  th*  Talue  of  T'  at  burnt  is 


"  am;  ''b 


.An.  . 

FK^W  • 


(10. 17) 


Ccnbining  equations  (10.16),  (10.17),  and  (lO.'S)  then  yields 

(1  +  j)  =  1  -  ^ _ ^J"|  f» 

[  2  A-i  Tj'J  ^  • 

If  a  near  parameter,  is  defined  as 


(10.10) 


s  y 


2{k-k^) 
“  ■  V  V  ■ 


(10.19) 


equation  (10.18)  finally  reduces  to 

(1  +  /)  ^*  =  1  -  T'  (10.20) 

idiieh  Tery  canreniently  has  the  same  form  as  equation  (4.6)  with  V  replacing; 
V,  and  replacing  'f.  S<q?arating  variables  and  integrating,  X  is  obtained  as 
a  function  of  T* 


X  = 


r  - 1 


- 1 


(10.21) 


which,  of  course,  has  the  same  form  as  equation  (4.7). 


Similarly,  if  a  new  pressure  variable,  P' ,  is  defined  as 


*  »  *  I 

the  relationship  between  />*,  F*,  and  is  found  to  be 

f  +  1 

P’  =  T »  =  T'  (l  -  (10.-23) 

which  has  the  same  fora  as  equation  (4.9).  The  value  of  F*  at  which  P' 
attains  its  mathematical  maximac,  P^ ,  is  thus 

V  .=  -i,  (10.S4) 
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and 


f '  +  1 


r-.  . 


f-.±l 


(10.2S) 


Tbe  value  of  T  is  found  by  substituting  the  new  dimensionless  parameters 
into  equation  (10.8) •  Thus, 


(10.26) 


which,  although  not  in  the  same  form  as  equation  (4. 16) ,  reduces  to  that 
equation  when  .=  y. 

A  special  case,  which  till  now  has  not  been  considered,  occurs  when 
•p'  .=  1.  Although  the  value  of  y  is  always  greater  than  one,  f  is  not  so 
restricted.  (We  shall,  however,  assume  -f'  >  0  since  otherwise  the  graphs 
become  quite  unwieldy.)  Letting  p'  .=  1,  equation  (10.20)  becomes 


dV  _  1 

'd'f  ■  r"i 

so  that  7*  and  X  are  related  by  the  expression 

¥•  =  in  (1  +  X) 

and,  P'  and  V*  by  the  expression 

p.  =  V' 

^  ''  dx  * 


(10.27) 


(10.20) 


(10.29) 


The  solution  is  now  complete.  Figiires  2  and  3  which  relate  the  vari¬ 
ables  X,  V,  P,  and  y  for  the  case  C  =  \L,  may  be  used  to  relate  the  new 
variables  X,  F',  P' ,  and  -p'  for  the  case  C  •=  \j  L  +  K2  -  Values  of  -p' 
need  not  be  greater  than  one,  hence  the  inclusion  of  values  of  y  <  *1  in 
the  figures. 

XI.  G31AIHICAL  DEIEEMENATICIN  OF 
MAXmW  FEESSUEE  AND  EJECTICK  VELOCITY;  C.=  \jZ,  +\,l2 


The  determination  of  maximum  pressure  and  ejection  velocity  for  the 
case  of  a  quadratic  form  function  is  considerably  more  involved  than  that 
for  a  linear  case.  This  is  not  surprising  since  "p' ,  unlike  the  constant  y. 
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is  a  fimction  of  the  maximum  pressure.  However,  in  the  case  of  extremely 
regressive  or  progressive  grain  configurations,  the  method  developed  here 
shoiild  lead  to  more  accurate  results  than  does  the  method  of  Section  VII. 

Bie  following  parameters  are  defined  similarly  to  those 

in  Section  VII: 

(11.1) 

^  {t\)  ‘'i 

(11.2) 

(11.3) 

Q'l  =  <7*,  /a^ 

(11.4) 

and 

O' 2  ~  7 *2  * 

(11.5) 

and 


The  dimensionless  variables  P\  F',  and  F'j,  are  then  written, 

P' 

V  =  Q\v  .=  q'  ^v/ <X 


(11.6) 

(11.7) 

(11.8) 


P'  t  P 

Since - .=  — -s-  ,= - ,  Figures  6a  and  6b  may  be  used  to  determine  a  as  a 

<?'i  «7i 

P' 

function  of  n  and  the  ratio  . 

7  1 

The  procedure  for  determining  maximum  pressure  and  ejection  velocity 
for  a  charge  with  a  quadratic  form  function  may  now  be  outlined. 
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4b 


Step  1 

a.  Given  the  firing  data  calculate  the  values  of  the  parameters: 
X  = 

Xj  =  X  -  p(Sfc  -  So)  /2 

-  Cfc/P 

r  =  (1  +  K^)f 

*c 

C*  -  ^ 

^  1  i2t  / 

.  A 
■  «?X" 

'  (t  )  Fxfir  » 

12iU„ 

-U--  • 


‘J'j  = 
0*3  = 


b. 

Assune  •f' 

=  y. 

c. 

Calculate 

II 

d. 

Determine 

=  / 

Step  2 

-a. 

Assume  avalue  of  P 

value  of  P*  from  Figure  3  corresponding  to  the  smaller  of  the  tvio  values 
as  determined  in  Step  1. ) 

b.  Calculate  • 

/  P'  \ 

c.  Determine  a  =  /  r,  *7^)  from  Figure  6. 

\  ^  V 

51 


53 


HWL  HSPOBT  m.  1758 


F  =  388,'250  ft  lby/lb„ 

Jj  =  0.089 
iTj  =  0.25 

=  0.06361  in 
n  =  0.372 

So  =  67.67  in*,  Sf,  .=  97-28  in*- 
Uc  =  159. 1  in* 
w  =  360  lb„ 
ar„  =  3.333  ft 
0  =  70.5® 
y  =  1.24 

p  =  0.05661  lb„/in* 

7)  =  30  in*/lb„ 

determine  the  maximum  pressure  and  ejection  velocity. 

The  solution  is  shoim  in  Table  5.  The  results,  -  1684  lb y-,  in* 

and  V*  =  56.8  ft/sec  compare  favorably  with  1647  Ibf/  ii»*  and  61.6  ft/sec 
obtained  from  numerical  integration.  (Disregarding  \2,  and  using  the  method 
of  Section  VII  yields  the  results  P^^^  =  1930  psi  and  =62.1  ft/sec.) 

XII.  rESIOf  OF  A  SCALE  UGOIEL 
TO  REPRDIUCE  BALLISnCS  OF  A  LARGER  CAD 

The  dimensionless  variables  defined  previously  may  be  xised  to  determine 
families  of  CAD's  having  id^itical  pressure-distance  and  velocity-distance 
ciurves  provided  the  distance  scales  are  suitably  adjusted.  A  very  useful 
application  of  this  feature  is  that  it  permits  the  determination  of  the 
dimensions  of  a  small  scale  model  to  reproduce  the  ballistics  of  a  larger 
parent  CAD.  The  use  of  a  small  model  as  a  testing  apparatus  could  lead  to 
considerable  saving  in  ei^>enditure. 
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all  such  devices  have  identical  ejection  velocities.  We  shall  call  those  CAD's 
satisfying  equations  (12.5)  through  (12.9)  "ballistically  8iniilar'\  If  the 
seme  propellant  ccnposition  is  used,  the  requirements  for  ballistic  similarity 
reduce  to 


=  VjT 

(12.10) 

h 

=  l/jT 

(12.11) 

Cb 

=  v^At 

(12.12) 

Uc 

-  v^At 

(12.13) 

W 

=  v^At 

(12. 14) 
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where,  again,  the  are  ocnstants.  Note  also  that  equationa  (12.  U)  and  (]£.  12) 
inply  that  the  foim  fimctioD  coefficients  of  ballistically  similar  devices  must 
satisfy 


Aj  =VgA 

(12.15) 

A]  =  v-i  A/r  ,  and 

(12.16) 

As  =  *^8  A/  t* 

(12.17) 

since  then  and  only  then  will  the  equation 

Cj  =  Lf,  +  Aj  Lj*  +  Aj  Lj,* 


be  valid. 

Ihe  choice  of  the  parameter  r  is  arbitrary  within  certain  limits.  It  can 
be  any  function  of  the  six  variables  A,  C^,  Uf.,  and  w  provided  the 

function  is  consistent  dimensionally.  For  practical  reasons,  however,  it 
should  be  a  simple  function  so  that  given  the  value  of  one  of  the  six  vari¬ 
ables,  the  values  of  the  other  five  may  be  determined  directly.  As  practical 
selections  of  r  we  might  hare: 

T  ~  A 


T 


(the  piston  dismeter) 


T  =  1/A  . 

The  following  numerical  exanple  will  help  to  illustrate  the  use  of  equa¬ 
tions  (12.10)  -  (12.14). 

A  CAL  has  the  following  dimensions: 

A  =  4.92  in* 

Uf.  =  159.1  in* 

w  =  300  lb,, 

X„  =  8.2  ft 

5  =  0»  . 
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Its  standard  chazige  is  eaif>osed  of  six  uninhibited  single-perf  grains  with 
dimensions 

D  =  0.51  in 
d  =  0.146  in 
h  -  z.  126  in 

idiich  lead  to  the  form  function  coefficients 
Aj  =  2.4541  lb„/in 
Aj  =  -1.4837  lb*/in2 
Aj  =  0  lb*/in3  . 

The  propellant  has  the  ballistic  properties: 

B  =  0.009436  (in/sec)  (in^/lby)" 

Cf,  =  0.21104  Ib^ 

Lf,  =  0.091  in 

n  =  0.56 
Ty  =  2410  •K 

p  =  0.06  Ib^/in* 

7,  =  30  inVlb„ 

y  =  1.239  . 

A  scale  model  is  to  be  made  with  a  cross-sectional  area.  A,  of  1  in^; 
vdiat  should  be  the  other  dimensions  of  the  device  and  those  of  the  propellant? 

We  shall  solve  the  problem  for  two  cases  by  specifying  two  different 
parameters  for  t  . 

Case  I:  t  =  A 

Substituting  the  value  r  =  A  into  equations  (12.10)  -  (12.17)  we  have 

Lb  =  vjil 

Cb  = 
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V  =  Vji* 

Ai  =  vgJ 
A]  =  V7 

A,  =  -vg  /i  . 

We  then  calculate  the  Talues  of  the  v  ^  from  the  given  dimensicna  of  tbe  parent 

CAD. 

*'l  =  8.2/4.9B  =  1.6667 

Vj  =  i^/ A. =  0.091/4.92  =  0.018496 
Vj  =Ci,tA*  =*0.21204/(4.92)*  =  0.0067184 
Vg  =  Oc/A*  =  159.1/(4.92)*  =  6.6726 
Vj  =  r/4*  =  300/(4.92)*  =  12.3934 
Vj  =  A| /A  =2.4541/4.92  =  0.49680 
«,  =  A,  =  -1.4837 
=  AjA  =  O  . 

Setting  A  =  1  f  ae  then  calculate  tbe  model  dimensions: 

=  VjA  =  (1.6667)  (1)  =  1.6667  ft 
Lf,  =  VjA  =  0.018496  in 
Cf,  =  V,  A*  =  0.0067184  Ib^ 

Oc  =  'VgA*  =  6.5726  in* 
w  =  V5  A*  =  12.^93  lb  ^ 

Aj  ='VgA  =  0.49880  lb  „/ in 
Aj  =-V7  = -1.4837  lb ^/in* 

A3  ~'V^/A  =0  1b^/in*  . 
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If  six  single-perf  gralsB  are  desired,  tbe  foim  function  coefficients  lead  to 
the  grain  dimenaions 

D  -  0.  365  in 

d  =  0.291  in 

h  =  0.635  in  .  * 

fbr  purposes  of  ccB^parison  the  dimensions  of  the  tao  devices  are  listed  in 
Table  6. 


Table  6  OdPARESCN  OF  mUENSIGNS  OF  PARENT  CAD 

AND  HEEL  CAD  OSING  t  =  A,  THE  PlgTCW  CROSS- 3E1CT1CKAL  AREA 


Parent  CAD 

Model  CAD 

A 

4.92  in* 

1  in* 

Cb 

0.21104'  Ib^ 

0.00672  lb 

h 

0.091  in 

0.0186  in 

fJc 

159.1  in» 

6.57  in* 

w 

300  Ib^ 

12.39  Ib^ 

% 

8.2  ft 

1.6667  ft 

M 

6  grains 

6  grains 

D 

0.51  in 

0.366  in 

d 

0.146  in 

0.291  in 

h 

3. 126  in 

0.636  in 

*Tbe  grain  dimensions  are  calculated  from  tbe  fomulas  given  in  i^pendix  A, 
i.  e. , 

Cb  =  (0*  -  * 

^2  = 


for  /f  =  6. 
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Usiog  nunerieal  integx*actioD,  tha  ballistic  perfoznance  of  the  two  devices 
was  calculated.  Thair  pressure  versus  x/r  and  velocity  versus  x/r  curves 
were  found  to  be  identical,  m  particular,  they  teth  yield  Pm^-r  ~  1403  psi 
and  =  101.2  ft/sec.  The  results  are  plotted  in  Figure  9. 

Case  II:  t  =  fl  =  2  »  the  Piston  Dlasieter 

Substituting  the  value  r  .=  0  into  equations  (32.10)  -  (12.17)  we  have 
*«  =  I'l 
=  vj  n 
Cf,  =  v,n» 

Uc  = 

w  =^50* 

Aj  =  n* 

A  j  =  fl 

As  =  >^8  • 

Substituting  the  given  dinensionB  of  the  parent  CAD  into  tte  foregoing  equations 
we  then  calculate  the  values  of  the  . 

=  8.2/2.5029  =  3.2762 

Vj  =  ij/n  =  0.091/2.5029  =  0.03636 

V,  =C{,/n’  =0.21104/(2.5029)®  =0.01346 

=Uc/n*  =159.1/(2.^329)®  =0.01346 

Vj  =i.;/n®  =  300/(2.5029)®  =10.147 

Vg  =  A|  /fi2  =  2.4541/(2.5029)  ®  =  0.  3917 

V7  =  Aj/n  =  -1.4837/2.5029  =  -0.5928 

*^8  =  A3  =0 

Setting  n  =  1.1284  in,  for  which  A  =  1  in^  ,  the  model  dimensions  are  then 
calculated. 

x„  =  VjO  -  (3.2762)  (1.1284)  =  3.697  ft 


h  =^20  =  0.04105  in 
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Cfc  -  v, fi»  =  aoia34  ib^ 

Oc  -  =14.80  in* 

V  ='V^C^  =’27.49  lb„ 

Ai  =  ^^  0*  =  0.4908  lb  ^/In 

A,  =’v^n  =  -0.6609  lb*/ in* 

Aj  =  Vj  =  0  1b* /in* 

Taking  M  ■—  6  grains^,  tha  grain  dimenBions  are  calculated  m»<«g  the 
fomulas  on  page  2  of  Appendix  A,  to  be 

D  =  0.'230  in 
d  =  0.0680  in 
h  .=  1.410  in  . 

Ualng  mnerlcal  .Integration  tecbnlques,  the  ballistic  perfonnance  was  calcu¬ 
lated  for  the  model  device.  As  before,  the  parent  CAD  and  the  model  CAD  were 
found  to  have  identical  preasure  and  velocity  versus  x/r  curves.  The  conpara- 
tive  dimensions  are  listed  in  Table  7. 

.^As  in  the  pnaviobs  case,  M  may  be  set  at  any  value  for  idiicb  D,  d,  and  h, 
have  positive  values.  If  M  =  -3,  for  exen^le,  the  dimensions  are  calculated 
to  be 

D  =  0.370  in 
d  =  0.214  in 
h  =  1.  410  in  . 

The  Isportant  thing  to  reoenber  is  that  the  charge  configuration  is  not 
unique.  Ahy  configuration  that  has  the  form  function, 

C  =  0.4908  L  -  0.6689  L* 

and  web 

Lj,  .=  0.04103  in 


will  produce  the  sane  resxilts. 
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Table  7  < 

1 

Parent  CAD 

Model  CAD 

A 

4.92  in* 

1  in* 

Cb 

0.21104  lb. 

0.0193  lb. 

h 

0.091  in 

0.041  in 

h 

159.1  .in* 

14.56  in* 

SI 

300  1b, 

27.49  lb. 

^w 

8.2  ft 

3.70  ft 

M 

6 

6 

D 

0.51  in 

0.23  in 

d 

0.146  in 

0.066  in 

h 

3.126  in 

1.41  in 

Tran  the  results  of  Cases  ;i  and  II  ae  conclude  tiiat  it  .is  possible  to 
detenniiie  the  dimensions  of  a  small  scale  model  that  is  baUistioally  similar 
to  a  larger  parent  CAD.  The  model  could  be  \)sed  for  e;q>erimental  purposes. 
Once  the  proper  charge  aeigbt  and  aeb  size  were  determined  for  the  model,  the 
proper  charge  design  for  the  parent  device  coxild  be  deduced  from  the  results. 

If>  now,  -r  is  taken  to  be  x^,  all  devices  with  equal  ratios, 

^b/^c  ’  will  satisfy  equations  (12.10)  and  (12.13). 

Ibr  a  particular  device  the  remaining  equations  determine  the  charge  design 
and  accelerated  mass.  Is  can  then  say  that  for  a  given  CAD,  any  device  with 
the  same  ratio,  Off/Ogf  may  be  used  as  an  ezperimentail  model.  The  follow¬ 
ing  .nvmerical  example  .illustrates  the  principle. 

Device  I  has  the  following  dimensions: 

A  -  4.92  in* 

Uc  ~  159.1  in* 
w  =  300  1b^ 

-  8.2  ft 


e 


=  0*  . 
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Tte  propellant  haa  the  prcpertiea: 

B  =  0.009436  (in/see)  (ii^/lby)" 

Cfc  =  0.21104  Ib^ 
ifc  =  0.091  in 

n  =  0.58 
Ty  =  2410*  K 
p  =0.06  1b^/in* 

7}  =  30  in»/lb„ 
y  =  1.239 

and  fonn  function  coefficients: 

A|  =  2.4541  lb„/in 
Aj  =  -1.4837  Ib^in* 

Aj  =  0  lb„/in3  . 

The  resulting  pressure  and  velocity  versus  x/jt^  curves  are  shown  in  Figure  10. 
TWO  other  devices  have  dimensions: 

Device  II: 

A  =1  in^ 

Uf.  =  3.944  in3 

=  1  ft 
0=0* 

and 

Device  III; 

A  =  '2  in* 

=  23.661  in* 
x„  =  3  ft 

e  =  0*. 

(Note  that  Uij/Uf.  =  3.043  for  the  three  devices.)  Using  the  same  type  of  pro¬ 
pellant^,  what  grain  design  and  accelerated  mass  should  be  used  in  each  device 
so. that  all  three  devices  produce  the  curves  shown  in  Figure  10? 
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Beplaeing  T  by  .in  equatiODs  (22.10)  -  (12.17)  ns  solTe  for  the-v^ 
using  tbe  propertioe  of  Derioe  ;i 

■V,  =  Lj/j*;,  =  0.0m0 

•'t  =  =  0.008231 

•V5  =w/Ax^  =  7.4360 

Vt  =Ai/i  =0.4968 

=AjX,/i  =-2.4729 

*'•  “  A*  x**  /a  =  0  » 

froD  Mhieh  m  detemine  tbs  xalues: 


Device  :il 

Device  III 

h  ^ 

O.OmO  in 

0.03330  in 

Cb  = 

0.000231  Ib^ 

0.031386  lb. 

w  =VsAx^ 

7.4360  lb. 

44.6160  lb. 

« 

II 

< 

0.4988  lb, /in 

0.9976  lb,/ in 

A,  =  v^Alx^ 

-2.4729  lb,/in* 

-1.6406  lb,/ in* 

A,  =’V,A/x^* 

0  lb,/ in* 

0  lb,/in* 

A 

1  in* 

2  in* 

Vc 

3.944  in* 

23.661  in* 

1  ft 

3  ft 

Integrating  tbe  ballistic  equations  numerically  and  plotting  tbe  results 
versus  */x^  for  both  devices  produces  the  curves  sboim  in  Figure  10. 

At  this  point  seme  mention  of  an  apparent  relaticnship  between  the 
action  times  of  ballietically  similar  devices  should  be  made.  As  a  con¬ 
jecture  we  state  that  at  equal  values  of  x/r  ell  such  devices  have  equal 
valxies  of  t  /'T.  Ibe  conjecture  has  not  been  proved,  but  there  aie  strong 
axgixnents  in  its  fkvor.  Consider  tbe  following. 
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Tbs  previcnis  derlTations  indicate  that  erery  menber  of  a  family  of  ballia- 
tically  similar  devices  bas  the  sane  velocity  equation 


idiere  /  (x/r)  is  a  particular  function  of  x/r.  Separating  variables  the  equa¬ 
tion  may  be  integrated  to  yield 


0 


The  integral  is  improper  and  may  not  exist  since  /(o)  =  0.  Ignoring  the  fact, 
however,  And  assuning  the  existence  of  the  integral  the  stated  conjecture 
follows. 

To  substantiate  the  argument,  in  Figure  11  we  have  plotted  the  valiie  t/r 
veravis  x /t  for  Devices  I,  II,  and  III,  where  t  =  x,, .  (The  data  were  obtained 
nixnerically. )  Although  there  are  actually  slight  differences  between  the 
curves  they  are  \indetectable  from  the  figure  and  are  possibly  due  to  inherent 
errors  in  the  numerical  integration  method.  At  any  rate  the  evidence  seems  to 
indicate  that  the  conjecture  given  above  is  at  least  a  very  good  approximation 
to  the  real  case. 

The  conjecture  may  be  stated  in  another  form:  for  ballistically  similar 
devices  the  time  required  to  reach  the  point  x/r  is  proportional  to  r.  Stated 
in  this  fonn  it  is  easy  to  deduce  a  very  useful  application.  From  a  given 
preastire-time  trace  for,  say.  Device  we  can  detennine  the  corresponding 
trace  for  any  other  ballistically  similar  Device  Dj  simply  by  multiplying  the 
abscissa  values  by  the  ratio  •  The  same  is  likewise  true  for 

pressure-distance,  velocity-time,  and  velocity-distance  curves. 

In  our  discussion  of  ballistically  similar  devices  we  have  neglected  to 
consider  non-horizontal  firings.  This  has  not  been  unintentional.  A  study 
of  the  non-horizontal  equations  was  undertaken  but,  unfortunately,  no  single 
similarity  relationships  could  be  found.  However,  the  foregoing  analyses 
seem  to  be  approximately  applicable  to  the  non-horizontal  case.  Perfoimance 
calculations  were  made  for  Devices  I,  II,  and  III  assiming  6  -  90”.  The 
results  are  plotted  in  Figure  32. 

Although  the  velocity  curves  are  very  similar,  the  pressure  and  time 
curves  vary  by  as  much  as  Tfa.  This  would  indicate  that  in  the  case  of  non- 

horizcntal  firing  the  qaations  should  be  used  with  some  discretion. 
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TOHM  fUNCnOWS  OF  SCME  CCUIOM  GRAIN  CONFIGORmCWS 

Ttie  law  of  btiming,  known  as  Piobert's  law,  is  generally  adopted  by 
interior  ballisticians.  According  to  the  law  the  surface  of  each  grain  in 
the  charge  recedes  parallel  to  itself  as  buming  proceeds.  Some  ccmfinna- 
tioQ  that  the  law  .is  nearly  obeyed  is  obtained  by  firing  charges  8\ich  that 
the  propellant  is  not  all  burned  at  the  end  of  the  stroke.  The  shepe  of 
each  tjnhumed  grain  is,  in  general,  found  to  be  very  well  preserved.  In 
other  words,  simultaneous  ignition  of  the  surfaces  of  the  grains  and  buming 
by  parallel  layers  seem  to  be  reasonable  assumptions.  (See  Hunt  (e)  page  40. ) 

A  grain  whose  surface  increases  as  buming  proceeds  is  called  "pro¬ 
gressive”;  if  the  surface  decreases  it  is  called  "regressive";  and  if  the 
surface  remains  constant  it  is  called  "neutral".  Some  of  the  more  comnon 
grain  shapes  and  their  form  functions  are  described  here. 

1.  ^lindrlcal_Graini^ 

To  obtain  the  form  function  of  a  charge  made  up  of  hcmogenecus 
cylindrical  grains,  the  total  charge  weight,  Cj,,  is  equated  to  the  nuniber 
of  grains,  M,  times  the  density  of  the  propellant,  p,  times  the  volume  of 
one  grain.  Thus, 

Cb  =  Hp--f-h  (Al.l) 

where,  D  and  h  are,  respectively,  the  initial  diameter  and  the  initial 
length  of  the  grain.  After  buming  through  a  distance,  L,  the  charge 
left  iinbumed  is 

Cf,  -  C  =  M  p  ^  (D  -  2L)^  [h  -  ZL)  .  (kl.Z) 

Combining  the  two  equations  and  rearranging  terms,  the  form  function  is 
obtained  as 

C  =  +  'Zh)L  -  ’£  [ZD  ^  h)L^  +  41^  j  .  {A1.3) 

The  charge  is  regressive.  The  value  of  to  be  used  with  a  quadratic 
form  function  is 

a.  f = 
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:if  D  >  h,  tbe  web  is  giren  by 

ifc  =  -g-  .  (A1.5) 

If  D  '<  A,  tbe  web  .ie  given  by 

ifc  =  l|  * 

A  comnon  variant  of  the  grain  deaign  is  to  "inhibit"  the  cylindrical 
surface  so  that  burning  takes  place  only  on  the  circular  .ends  of  the  grains. 
A  neutral  charge  is  obtained  in  this  case,  the  fonn  function  being 

C  =  L  (A1.7) 

and  the  web  is  given  by  equation  {A1.5) 


2.  Single-Perf  Grains 

A  "single-perf"  grain  is  a  cylindrical  grain  pierced  by  a  cylindrical 
perforation  through  tbe  center  of  the  grain  parallel  to  tbe  length  of  tbe 
grain. 

Equating  the  charge  weight  to  the  nunher  of  grains,  jV,  tines,  tbe 
density  of  tbe  propellant  times  the  volune  of  one  grain  yields 

Cb  =  fl  (D*  -  <i*)A  (A2.1) 

where  D,  d,  and  A  are  tbe  initial  values  of  tbe  dieneter  of  the  grain, 
the  dianeter  of  the  perforation,  and  tbe  length  of  the  grain,  respectively. 
After  binning  through  a  distance,  L,  the  charge  left  unbumed  is 

Cf,  -  C  =  M  p  -  2L)  *  -  (d  +  2L)  » j  (A  -  .  (A2.'2) 

Combining  the  two  equations  we  have 


•2 
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Ibe  grain  .la  regreaslTe,  only  sli^itly  so  if  k  '>  >  D.  Since  the  foim  fXmction, 
equation  (A2.-3)  I  is  quadratic  the  ralue  of  X.|  .ia 


=  TiK^So  [*  ^  ] 


(A2.4) 


The  neb  .la 


or 


,  _  D  -  d 

h  -  ”4- 
I  -  * 


(A2.5) 


whicheTer  ia  the  amaller. 


If  the  ends  of  the  grains  are  inhibited  a  neutral  burning  charge  is 
obtained  with 


L 

h 


(A2.6) 


and 


r  .D-d 
h  -  ”4 - 


(A2.7) 


If  all  the  surfaces  except  the  perforation  are  inhibited  a  progressive 
charge  is  obtained  with  the  form  function 


c  =  (di  +  i2) 

/)*  -  d* 


(A2.8) 


with  X|  being  given  by 


X  - 


(A2.9) 


and  web, 


.  .D-d 

h  —~2' 


(P2. 10) 


3.  Seven-Perf  Grains 


A  "seven-perf”  grain  consists  of  a  fairly  short  cylinder,  pierced  by 
seven  equally  spaced  cylindrical  perforations  parallel  to  the  length  of  the 
cylinder.  It  is  generally  specified  that  the  diemeter,  d,  of  each  perfora¬ 
tion  is  approximately  one-tenth  the  diameter,  D,  of  the  grain.  The  appear¬ 
ance  of  the  end  section  is  shown  in  Figure  13a. 
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Tte  .initial  distance  betiieen  any  two  perforations  and  between  ai^  of  the 

D  -  ’Sd 

outer  six  perforations  and  the  curved  surface  of  tbs  cylinder  .is  — ^ —  .  :it 

D  —  *3d 

will  be  noted  that  when,  a  distance  been  burned  throu^  there  are 

o 

twelve  curvilinear,  triangular  prians  called  "slivers"  remaining  unbumed. 
(See  Figure  13b.)  Strictly  speakii^,  then,  a  seven-perf  grain  has  two  form 

D  —  'Sd 

functions-one  valid  for  0  <  L  <  -~s —  ,  the  other  valid  during  burning  of 

the  slivers.  In  practice,  however,  since  the  grain  is  about  85  per  cent  con- 
suned  at  the  moment  of  slivering  and  since  the  latter  fonn  function  is  unduly 
coofjlicated,  it  is  generally  assuned  that  the  initial  fonn  function  is  appli¬ 
cable  during  the  entire  burning  regime. 

Equating  the  total  charge  weight,  C^,  to  the  mmber  of  grains,  F,  times 
the  density  of  the  propellant,  p,  times  the  volume  of  one  grain  gives 


Cb  -  H  (P*  -  7d»)A 


(A3.1) 


irtiere,  as  before,  h  is  the  length  of  the  grain.  After  burning  throu^  a 
distance,  L$  the  charge  left  unbumed  is 


Cf,-C  =  M  p 


t[«>- 


S£)*  -  7(d  +  2i)* 


2L) 


(A3.2) 


Combining  the  above  equations  the  form  function  becomes 

C  = 


[h[D*7d)+-^  (P*  -  7d*)]i  +  [6A  -  2(P  +7d))L*  -  12L* 


I  .(A3. 3) 


=  — j| 

h{D^  -  7dJy  Y 

The  charge  is  progressive.  The  value  of  X|  for  use  in  a  quadratic  form 
function  is 

i;  -  ^b\j,Tzya-t  zh  A“(ir*‘-7d*)J  ■ 


TO  be  consistent,  it  becosies  necessary  to  define  the  web,  as  the  real 
root  of  equation  (A3. 3)  when  C  =  C^.  In  the  general  case,  for  which  A  >  2L^, 
we  then  have 

^  P-/7d  ^  P-2.6456d 

■"7.2916"  (^-5) 


i,  =  - -I 


•2(1-H^) 


otherwise 


=  i- 
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(A3.6) 
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bOBHBgnCW  FACIOK  TOR  RB3I3E/glCg  DDB  TO  SFRIWS 

Occasionally,  as  in  seme  "stores  separation  devices",  tbe  resistance  of 
a  spring  must  be  overcome  in  addition  to  tbe  other  resisting  forces  actii^  on 
the  accelerated  mass.  In  most  cases,  this  resistance  .is  quite  small  and  can 
be  accounted  for  by  a  slight  change  in  the  potential  eneigy  correction  factors 
derived  in  Section  VI. 


Let  I  be  the  known  spring  constant.  The  equation  of  motion  is 


,=  »  ,  si.  aj 


I  --f-  +  f  sin  0 1  +  /* 
and  the  energy  balance  equation  is 


(B.1) 


CyC{Ty  -  T)  =  "I  sin  e)  +  .  (B.2) 


If,  now,  tbe  tern  V  g  sin  6  +  lx  ia  taken  as  being  proportional  to  we 
obtain  the  same  fonn  for  tbe  co 

obtained  in  Section  VI,  namely 


obtain  the  same  fonn  for  tbe  correction  factors  for  f  sod  v_  as  was 

itwrr  "« 


1  +  — -JL-  )a 
3  -  2n  rxDsx 


(B.3) 


and 


=  ll- 


1  -  n 


3  -  Zn  •'m  • 

The  value  of  f  j  is  slightly  changed,  however,  and  is  obtained  as  follows 


(B.4) 


From  tbe  assvmption  of  proportionality  we  have 
X 

r 

xAp 


J  (1  g  Bin  0  +  IxJdx  P 


T  dx  . 


(B.5) 


Substituting  p  ~  W  into  the  foregoing  equation  and  inte- 

A  dt^ 

grating  we  obtain 


W  gx  sin  9 


(B.6) 
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The  equation  must  be  satisfied  for  v  .=  and  x  =  ;  hence,  we  have  as  the 

value  of  JTj 


^3 


{B.7) 
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and 


X  =  (l2A/ai)x 
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The  distance  ▼ariable  JT  is  related  to  the  Telocity  variable  F',  by  the 
expression 


/  =  (1-  /iF*)'#*  -  1 


where 


•  so  that  equation  (c.  1) 


(C.4) 

(C.5) 


may  be  transfomsd  into  an  equation  relating  the  time  rate  of  change  of 
acceleration,  a,  to  the  variable  F*.  Thus, 


2  ^  ^  ^  ^ 

(1  +m)k'2  (1  -  ^Jy•)'  *'*1  . 


(C.6) 


Differentiating  the  above  equation  with  respect  to  F',  equating  the 
resulting  expression  to  zero,  and  solving  for  F'  leads  to  values  of  F'  for 
which  a  is  a  singular  point.  The  value  of  F'  for  which  a  is  a  maximum  is 
subsequently  found  to  be 


F.  =  1 

f '  L  Var'  +  1  J 


This  value  of  F'  leads  to  the  expression  for  maximum  a 

\3 


max 


— zrn—*—  'P\y  )  f  s  /  sec. 
Fc 


(C.7) 


(c.e) 


where  <p{y')  is  a  function  of  y'  only. 

Ass\]m 
expression 


As8\jming  that  is  a  mathematical  maximum  we  may  substitute  the 


^  Vr'A  2r'  /  12F\a5R*l/ 

into  equation  (C.8)  to  yield 


max 


(C.9) 


1"“  Fc  V 


(C.  10) 


Denoting  by  co{y')  and  adding  the  as  yet  unspecified  correc¬ 

tion  factor  k  then  leads  to  the  approximation  given  above. 
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The  approximation  for  y*  was  obtained  in  the  folloiring  manner,  y*  is 
related  to  y  by  the  expression 

r  =  .  (C.U) 

Combining  the  above  equation  with  equation  (C.9)  yields 


r' 


X\F  ^ 


(C.12) 


In  practice  y'  never  differs  from  y  by  more  than  a  few  percentage  points  so 
that  we  may  approximate  the  expression  on  the  right  by 

(C.13) 

\ 

This  then  leads  to  the  y*  approximation 


(C.14) 


A  comparison  of  values  of  a  calculated  from  the  above  expression 

4 

with  those  obtained  experimentally  indicates  that  the  value  k  =  —  for 

both  the  Mk  327  Test  Set  and  the  Dahlgren  ballistic  pendulum  gives  gener¬ 
ally  good  agreement.  The  results  of  the  comparison  are  tabulated  below. 


Experimental  o  Calculated  a 

max  max 


Mk  327 

260.1 

g's  /  sec. 

254.9 

g's  /  sec. 

Test 

284.4 

g's  /  sec. 

313.3 

g's  /  sec. 

Set 

222.7 

g's  /  sec. 

241.4 

g's  /  sec. 

Dahlgren 

19.4 

g's  /  sec. 

18.9 

g's  /  sec. 

Ballistic 

13.0 

g'B  /  sec. 

12.8 

g's  /  sec. 

Pendulum 

16.5 

g's  /  sec. 

16.0 

g's  /  sec. 

A  "better"  value  of  k  might  be  obtainable  from  more  extensive 
enpirical  data. 
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APPUCmON  OF  IHE  LAEKaR  ROOT  OF  EgJATICN  8. 14 

For  one  reason  or  another,  the  ratio  of  barrel  volvine,  U^,  to  chamber 
volune,  Uf.,  for  most  CAD’s  is  of  the  order  of  1.  For  exanple,  the  ratios  for 
the  Mk  327  Test  Set  and  the  Dahl^ren  ballistic  pendulum  are  1.24  and  1.57 
respectively.  For  ratios  of  this  order  of  magnitude  the  methods  of  Sec¬ 
tion  VIII  are  tisually  sufficient  to  determine  charge  vveight  and  web  size  to 
meet  given  perfonnance  requirements.  Recently,  however,  a  device  was 
encountered  with  a  ratio  Uf,/  Uf.  =  5.  In  attempting  to  obtain  a  grain  design 
for  the  device  it  was  found  that  no  feasible  solution  could  be  obtained  for 
any  value  of  E  through  the  use  of  Figure  7a.  An  investigation  was  under¬ 
taken  to  determine  the  reason  for  the  failure  and  the  means  of  correcting 
the  situation.  The  results  of  the  investigation  follow. 

In  Section  VIII  B  it  is  shown  that  for  a  given  value  of  £  a  feasible 
solution  exists  if  and  only  if  satisfies  the  inequality 


The  derivation,  not  given  in  Section  VIII  B,  is  as  follows: 

^  Vb 


\0iBF\)  -  [oB)  FXW 

Vvl  <  FKLhV„i  =  FCf,V„i  =  Fp  -  i/^)  V, 


ml 


1  p^ij  It  _ 

2  -  2  ‘'P^c  n 
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lAere 


F„l  is  a  root  of  the  equation  * 

t'-5  =(i- -!• 


(D.2) 


(D.3) 


and 


E  = 


1  , 
^max^'a 


(D.4) 


Recall  that  equation  (D.2)  baa  two  real  roota,  the  amaller  of  irtiich  ia 
denoted  by  .  We  ahall  denote  the  larger  by  F,,l  . 


Tbe  bracketed  term  in  relation  (D.1)  must  be  poaitiTe,  bence  a 
neceaaary  condition  fbr  tbe  exiatence  of  a  feaaible  aolution  ia 


<  /,  .  (D.5) 

c 

Referring  to  Figure  7a,  tbe  value  of  F^j  corroaponding  to  is  found 

5 

to  be  approximately  equal  to  -  y  .  Subatitution  of  thia  v^ue  into  equa¬ 
tion  (D.3)  with  subaequent  use  of  tbe  relationabipa  (D.5)  and  (D.  l)  then 
yielda  the  result: 


If 


Oc 


>  R 


max 


idiere 


oil-  if 


(D.6) 


(D.7) 


{D.8) 


2 


the  larger  root,  ,  of  equation  (D.2)  must  be  used.  If  neither  (d.6)  nor 
(D.8)  is  true  the  analler  root,  ,  should  be  used. 

In  light  of  the  abore,  a  slight  modification  in  the  methods  of  Sec¬ 
tion  VIII  is  required.  In  Figure  7b  ne  have  plotted  versus  E  and  y . 
Notice  that  the  right-hand  menber  of  the  inequality  (D. l)  is  a  monotonically 
decreasing  function  of  E  when  F,,|  is  used,  while  it  is  monotonically  increas¬ 
ing  frtien  is  employed.  The  modifications  to  be  made  are  then  as  follows: 

In  Step  2  of  Parts  A  and  B  of  Secticm  VIII  first  detennine  once  and  for 
all  from  the  relationships  (D.6)  and  (D.8)  which  root  should  be  enployed  — — 

similarly  in  Part  C  except  that  ^  V  in  (D.8)  is  replaced  with 
The  rest  of  Part  A  remains  unchanged.  In  Step  2c  of  Part  B  proceed  as 

W 

directed  unless  F^j  is  used.  In  that  case,  if  ^  °  smaller 

value  of  E  and  repeat  Step  2.  Similarly,  in  Step  2d  of  Part  C,  if  E^  > 
select  a  smaller  value  of  E  when  using  F„,|  . 

The  above  procedixre  was  used  in  the  case  of  the  previously  mentioned 
device  where  Uf,/  Ug  -  b.  A  feasible  solution  was  found  using  and  veri¬ 
fied  by  the  numerical  integration  of  the  ballistic  equations. 

An  interesting  observation  may  be  made  from  the  relationship  (D. l). 

If  the  inequality  is  replaced  with  equality,  burnt  occurs  at  the  moment  of 
con^letion  of  the  stroke,  since  then  Fj,  =  F„j  .  For  a  given  propellant, 

barrel  volume,  ejection  velocity  and  maximum  pressure,  this  fact  then  enables 
us  to  detennine  two  charber  volumes  and  two  constant-surface  charge  configu¬ 
rations  to  give  the  desired  performance.  Furthermore,  burnout  in  both 
instances  should  occur  at  end  of  stroke.  Consider  the  following  example. 

The  propellant  has  the  following  ballistic  properties: 

B  =  0.01  (in/sec)'  (in*/lby)  ” 

n  =  0.5 

F  =  279,560  ft*lby/lb„, 

■y  =  1.2 
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p  .=  0.06  Ib^/in* 

7).=  16.67  in*/ lb*  • 

The  device  has  dimensions: 

=  40  in* 

^,  =  JTj  =0 

M  -  3106.1  slugs 

=  8  ft;  =  12  il  =  3840  in* 

0  =  0°, 

Determine  two  chamber  volumes  and  their  corresponding  constant-surface  charge 

configm*ations  to  yield  P  .=  1354  psi  andv_  =  15  ft /sec.  Burnt  should 

max 

occur  at  ejection. 

Solution 


Using  equality  in  the  relation  (l)  we  can  solve  for  Ug  to  obtain 


= 


Since 


Oi 


equation  (D.9)  in^lies  that 

^  •'‘'in* 


(D.9) 


(D.  10) 


(D.ll) 


From  the  given  perfonnance  requirements  we  calculate 


=  0.8068 


(D.12) 


and  determine  a  •=  0.036  from  Figure  lb.  Using  Figure  7a  we  obtain  .=  0.9 
and  the  corresponding  value  .=  1.55  frcm  Figure  2b.  The  required  charge 


4' 


nel^t  is  foimd  to  be 


C5  =  “  2.779  Ibj,  ;  the  initial  free  volune  is 


^  =  2477  in*  , 


80 


that  iL  =  -^  +  Ui  -  2523  in* . 

V  O  • 


Buxnt  ooeiirs  at  end  of  stroke,  hence  equation  (4.1)  msy  be  used  to  detennine 
the  web, 

ifc  =  ™  V,  =  0.4196  in  . 

In  smnary,  then,  the  Tallies 

Uc  =  2523  in* 

Cfc  =  2.779  Ib^ 

4  =  0.4196  in 

should  produce  the  desired  results. 

Similarly,  using  Figures  7b  and  2b  we  obtain  =  1.8  and  X„  =  6.4. 
Bqploying  the  same  procedure  as  above  we  then  have 

Uc  =  623.2  in* 

B:  I  Cf,  -  1.390  lb„ 

Lf,  =  0.4196  in 

as  the  second  set  of  values. 


Using  numerical  integration  techniques  the  ballistic  performance  was 
calculated  for  the  two  sets  of  values,  A  and  B.  The  results  are  shown  in 
Figure  15.  Notice  the  radical  difference  in  the  shape  of  the  curves,  though 
they  produce  nearly  the  same  p  and  v_.  The  distance  burned  through  the 

nicuv  ™ 

web  in  case  A  is  found  to  be  .3978  inch,  while  that  in  case  B  is  .3586  inch. 
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Hie  above  process  could  be  repeated  to  determine  the  required  chamber 
volune  for  p  fixed  at  1354  pai  and  various  values  of  v_ .  By  allowing 

to  take  on  values  between  14:8  ft/ sec  and  15.6  ft/sec,  the  conplete  range 
of  E  given  in  Figures  7a  and  7b  may  be  covered.  Numerical  integration  tech¬ 
niques  would  produce  a  large  family  of  cvurves,  two  examples  of  vdiich  are 
shown  in  Figure  15.  Depending  on  the  desired  pressure- time  characteristics 
an  "optionin"  chamber  volume  could  then  be  picked  out  from  the  family  of 
curves.  For  exanqile,  if  maxiimin  ejection  velocity  is  taken  to  be  the  cri¬ 
terion  for  "optiitum"  the  procedure  yields  an  optiimm  chandler  volume  of  about 

2710  in®  for  which  v_  =  16.4  ft/ sec  and  p  -  1375  psi. 

IB  max 

The  foregoing  suggests  the  possibility  that  for  a  given  device  with 
fixed  dimensions,  improved  performance  could  be  obtained  sinply  by  loading 
the  chamber  with  inert  material  and  decreasing  the  charge  weight. 


AH^EM)IX  E 
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REFINQfEllTS  TO  GRAIN  DESIGN  TO  RQIDER 
■«ND  V-  WITH  BDRNDOT  NEAR  END  OF  STROKE 

flWT  NS 

Several  methods  of  determining  constant-surface  charge  requirements  to 
meet  given  CAD  perfonnance  epee  if  icat  ions  are  discussed  in  Section  VIII.  One 
of  the  shortcomings  contnon  to  all  these  methods  is  the  fact  that  quite  often 
they  lead  to  designs  that  are  only  fractionally  consumed  at  the  end  of  the 
piston  stroke.  A  method  of  correcting  the  situation  is  given  in  Section  IX, 
but  the  correction  is  not  alvrays  practical,  especially  for  those  CAD's  having 
high  loading  densities.  Furthermore,  although  constant-surface  grain  designs 
are  usually  adequate  for  those  devices  having  barrel  volume  to  chsrnber  volume 
ratios  tas  1,  they  generally  perform  poorly  in  those  devices  with  ratios  >  >  1. 
For  these  reasons  we  retxim  to  the  problem  of  grain  design. 

Let  conditions  a)  and  b)  of  Section  VIII  be  replaced  by  the  stronger 
conditions 


a')  "Bumf  shall  occur  at  end  of  stroke,  i.  e. ,  Fj,'  =  Vf„\  . 

b')  p  shall  be  a  mathematical  maximum,  i.  e. ,  VJ,  >  VL  =  l/y' . 

Under  these  conditions  we  shall  attempt  to  deteimine  the  charge  weight  and 
grain  geometry  required  to  produce  a  given  maximum  pressure  and  ejection 
velocity. 

Condition  a')  permits  the  inmediate  determination  of  the  web-size  from 
equation  (lO. l) .  Thus, 


(E.  1) 


Enploying  the  definition  of  initial  free  volume  together  with  equation  (lO.  14) 
yields  an  equation  for  the  charge  weight, 

where  ~  ISAx^,  the  barrel  volume.  Condition  a')  also  permits  the  rearrange¬ 
ment  of  equation  (10.15)  to  yield  an  expression  for  ,  the  first  form- function 
coefficient. 
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UHL  Rggg’ 


^ '  F-ir^s  ■  ‘*-®’ 


Finally,  equation  (10.2)  nuy  be  used  to  detennine  the  value  of  Xj  ,  the  second 
coefficimit. 


X-Xj  Cff/L/f-K^ 

* "  '^r "  "”6””  • 


(E.4) 


The  requii^ed  charge  wei^t  and  grain  geometry  are  given  by  these  four  equations. 
Hoirever,  before  they  can  be  evaluated  we  must  first  detennine  F^|  ,  and  from 
the  given  firing  specifications. 


Gy  condition  a*,  is  related  to  and  y' ,  as  before,  by  the  expres¬ 


sion 


fi- 

f '  - 1  L 


y'  •  I 


(1  + 


]  .  r  1 


Vf  = 


In  (1  +  Jj  ,  y'  =  1  . 


(B.5) 


Using  a  procedure  identical  to  that  employed  in  Section  VIII  we  define  the 
piezometric  efficiency  by  ^ 

(E.6) 


and  the  parameter  F'  by 


F' 


y'  ♦  1 

— 


f  .  r  =  1 


X^,  Fp|,  and  y*  are  then  further  related  by  the  equation 


X  =  C-  F'* 


(E.7) 


^Nbte  that  from  equations  (10.25)  and  (10.29)  we  have  the  identity  relation¬ 
ship  F'  -  1/2  P’/f. 


2 


Tbus  far  wa  have  two  equations  in  the  three  unknoims  F^l  and  y'.  He 
require  one  acre.  If  we  define  a  "volunetric-efficiency"  as 

-  Vp%-  te-8) 

a  rather  lengthy  algebraic  manipulation  yields  the  third  equation. 


.  ji-fBini 

"  ■  ij-X'Z'  * 

i-^^i  uiij  ^ 


'The  derivation  is  as  follows: 


=  n 


(..L..y  „1  .(ty  .it.  r 

[aBFKJ  [ob)  f  \j  it 


■  ^1  Ki  FL(, 


=  X  ymi  FC^ 


F  p-U-  \  "2  -  Uyx-)  • 


Also,  from  equation  (10.19)  we  have 


21^  -  1) 


Solving  for  r-  yields 


x:  ‘  -I-  ^ 


Ihe  condiination  of  equations  (a)  and  (b)  then  leads  to  equation  (e.9)< 
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If  the  given  specifications  satisfy  ceirtain  restrictions  the  three  non¬ 
linear  equations  (E.5)»  (E.7),  and  (e.9)«  may  be  solved  for  the  variables 
1  and  y'  by  any  one  of  several  possible  numerical  techniques.  The 

equations  (E.I)  - (E*4)  may  then  be  evaluated  to  determine  the  grain  design 

requirements. 

The  mxnerical  solution  of  a  system  of  non-linear  equations,  however, 
can  be  quite  tedious.  Moreover,  there  is  no  guarantee  that  equations  (E.5), 

(E. 7),  and  (E.9)  possess  a  real  solution  for  a  given  value  of  E.  Conse¬ 
quently,  a  less  direct  but  sinpler  method  of  attack  is  suggested.  Before 
describing  the  procedure,  however,  we  shall  have  need  of  some  further  relation¬ 
ships  between  the  variables. 


In  order  for  a  grain  design  to  be  meaningful,  the  burning  surface,  S, 
must  of  course  always  be  non-negative.  In  particular  we  must  have 


So  >  O 

Sb>0 


(E.  10) 


where  and  ere,  respectively,  the  initial  and  final  burning  surfaces. 
From  the  form  function  equation  we  obtain 


_  dC  _  2. 

p  dl  p 


(X,  +2X2  I.) 


(E-11) 


so  that  the  required  inequality  relationships  become 

X,  >  0 

X|  2X2  ^  0  . 


(E.I2) 


These  inequalities  restrict  the  permissible  range  of  values  of  y'  and 


The  first  inequality,  Xj  >  0,  implies  by  definition  that  V'  must  be  non¬ 
negative  for  all  values  of  1.  F^j  must,  therefore,  be  positive.  Also, 
since  y'  =  l/v'ff,  y*  must  be  non-negative. 


=  0  implies  Xj  =0.  The  equations  of  Section  X  are  not  applicable  in 
this  case  since  F'  and  P'  are  then  undefined. 
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Hie  second  Inequality,  +  Z\.^  Lf,  >  0,  implies  that  the  ratio  X/X| 
must  satisfy 


A-  >  i  . 

X,  -  2  » 


(E.13) 


this  foUoHS  fron  equation  (E.4).  ]9Bi>loying  equation  (b)  in  footnote -2,  wa 
then  bare  another  restriction  on  y‘  and  i.  e. , 


^mi  (y'  -y)  ^  1* 


(E.14) 


Combining  the  above  restrictions  with  condition  b')  leads  to  the  follow¬ 
ing  conclusion:  we  need  only  consider  those  solutions  to  equations  (E.5), 

(E. 7) ,  and  (e. 9)  for  which 


y  >0 


— i-;  <  -z-  <  ¥•. 


---  <  VI,  <  — ^ — 


for  y*  <  y 


for  y  >  y  , 


(K.15) 


F^l  and  7'  are  said  to  be  in  the  "solution  domain"  if  they  satisfy  these 
inequalities  as  well  as  equation  (E. 15)> 

He  may  now  describe  a  suggested  procedure  for  determining  change  weight 
and  grain  geometry  to  yield  a  given  performance. 

Steg_l. 

Given  the  maximum  pressure,  p  and  ejection  velocity,  v_,  together 

0Q83C  ^ 

with  the  other  pertinent  ballistic  parameters  calculate  the  values 


^3  =  sin  0/v^ 

^max  "  ^max  ^  ^  3-^n  ) 

"  z~-Zn  ^») 
V  -  (1+^1)  w/yc 
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+  1 


E.= 


Ii'Uml.. 

^nax^** 


1  "  T 77/7 


fJh  = 


^1  = 

^2  = 


Ub 

U-'^eJUc 

_./.L_-. 

l-7£i 


Steg_^. 

a.  ABsixne  a  value  of  X^.  (initially,  take  =  Mj  •) 

b.  From  Figure  2b,  select  a  value  of  y'  and  F^j  in  the  solution  doraain 
corresponding  to  * 

c.  Calculate 


^The  solution  domain  is  determined  in  the  following  manner.  In  Figure  2b,  con¬ 
nect  those  points  for  irtiich  F*  =  ---  .  Also  connect  those  for  which 

y' 

F'  =  — —  •  The  solution  domain  lies  within  the  region  bounded  by  these 

y'  -  y 

two  curves.  In  Figure  16,  we  have  plotted  the  solution  domain  for  y  ~  1.295. 
He  shall  use  this  figure  in  place  of  Figure  2b  in  a  sait^le  problem  given 
later. 
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d.  If  the  aseaisd  value  of  eq^iale  the  ealeulated  value,  go  on  to 
part  e;  if  not,  repeat  parte  a.  —  e.  using  the  calculated  value  of  /^,  and 
the  sane  value  of  y*. 


e.  nroD  Figure  '3,  deteimine  the  value  of  corresponding  to  the  value 


of  y. 


f.  Calculate 


^  PViI  • 


Sten -3. 


a.  Repeat  St^  2  for  several  different  selected  valws  of  and  so 


as  to  cover  the  eolutlon  donain. 


b.  Plot  the  values  as  ootained  in  Step  *2,  E  and  versus 

0.  Select  the  value  of  and  F^i  corresponding  to  the  desired  valiw 
of  £  as  obtained  in  Step  1.  (Note  that  there  msy  be  one,  two,  or  no  solu¬ 
tions.) 


Step  4; 


a.  Determine  a  =  /  (n,  from  Figure  1. 


b.  Calcxilate 


h  -  -J- 


Cb-pPc 


V  - 


S  =  --1  ,  0  <  L  <Lb  . 
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These  equations  detemine  the  required  charge  weight,  web-size,  and  surface 
versus  buming-distsDce  relationship  to  yield  the  given  perfonnance. 

Sanjple  Problan  i 

Given  a  CAD  and  propellant  with  the  following  properties: 

A  =  4:92  in* 

B  =0.01919  (in*/lby)"  in/sec 
F  =  328,250  ft.  Iby/ lb  „ 

/j  =  0.089 
A'j  =  0.25 
n  =0.372 
Uc  =  1S9.1  in* 
u)  =  360  Ib^ 

=  13.474  ft 

e  =0* 

y  =  1.24 

p  =  0.05661  Ib^/ in* 

y\  =  17.665  in* /lb „  . 

Deteraiine  a  grain  design  to  yield  an  ejection  velocity  of  120  ft/sec. 
The  maxiinun  pressure  should  be  no  greater  than  1600  psi.  Burnt  should  occur 
at  or  near  ejection. 

Solution 

Proceeding  according  to  the  instructions  in  Step  1  we  calculate 
ATj  =0 

^max  ^  1600  psi 

Vjj  =  120  ft  /  sec 
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¥ 

12.185  slugs 

y 

= 

1.295 

E 

> 

0.827 

= 

0.02967 

Ob 

= 

796.5  in> 

= 

5.200 

= 

0.000856  . 

Oontinuing  on  to  Step  >2,  we  guess  .=  /Xj  .=  5. ’200.  Flxxi  Figure  16  we 

select  the  maximm  values  of  y‘  end  in  the  solution  donain  correspond¬ 

ing  to  =  5.2.  Thus,  -y*  =  0.45;  F^j  .=  2.35.  We  then  calculate  the 
value  of  frcm  the  equation 


,  _  5.2 

- - — 


1  -  0.03006 


7"2" 

^2.36 


and  obtain  5.265.  Repeating  the  process,  this  thne  guessing  X^  ~  5.265, 
there  is  no  perceptible  change  in  F^l,  so  that  the  values 

Vl  =2.36 
f'  =  0.45 
=  5.265 

are  taken  as  satisfying  equations  (E.5),  (E. 7),  dnd  (E.9)  for  some  value  of 
E-  From  Figure  >2  we  obtain  Pjf  .=  0. 64’  corresponding  to  y'  .=  0. 45  and  calcu¬ 
late 

If ',2 

E  =  -prf-  =  0.82. 

In  Step  3  we  repeat  the  process  for  different  values  of  y'  and  obtain  the 
following  table. 
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Table  8  SOLDTICKS  CF  EffJATICNS  (E.5), 

(E.7)  .  AND  (E.9)  FOR  VARIODS  VALDES  OF  E 


5.2  5.265 

5.2  5.249 

5.2  5.231 

5.2  5.218 

y‘ 

0.45  0.45 

0.6  0.6 

0.8  0.8 

1.0  1.0 

2.35  2.35 

2.22  2.22 

2.01  2.01 

1.8  1.8 

5.265  5.265 

5.249  5.249 

5.231  5.231 

5.218  5.218 

p'm 

0.64 

0.525 

0.43 

0.37 

E 

0.82 

0.89 

0.90 

0.84 

Since  we  are  only  interested  in  values  of  £  >  0.83  we  have  sufficient 
data.  The  values  are  then  plotted: 


1.0 


0.9 


0.8 


'mi 


E 


j^parently,  any  value  of  between  approximately  1.8  and  2.3  should 
meet  the  given  specifications.  Depending  upon  the  psrticxilar  choise  of 
the  corresponding  value  of  may  then  be  picked  off  from  the  graph.  For 
our  purposes  we  shall  select 

=2 

=  5.23 
E  =  0.9 


and  go  on  to  Step  4. 


K) 


Before  we  can  detemlne  a  we  calculate  the  value  of  p  from 

p  =  IliwI 

'^max  ■  E  Ax„ 

to  obtain  P  =  1470  pai.  From  Figure  6a  the  value  of  a  .=  0. 015  is 

IMJC 

obtained  for  n  =  0.372  and  p^^  -  1470.  Finally,  we  calculate  the  values 
^6  =  "w  =  0.06554  in 

-6  "  (l  -  7/:\)  =  0.3963  lb„ 

Vv  ^ 

^1  =  =  3.3246  lb„/in 

Xj  =  =  17.633  lb„/in2 

X  2X  L 

and  S  =  =  55.20  +  623.0  L  in*  0  <  L  <  0.06554. 

P  P 

These  values  provide  a  preliminary  design  which  may  be  modified  on  the  basis 
of  the  results  of  a  nvmerical  integration  of  the  ballistics  equations  and/or 
experimental  firing  results.  For  example,  using  numerical  integration  the 
above  design  is  found  to  yield 

^max  ~ 

=  112  ft  /  sec  . 

In  addition,  the  web  is  fcvmd  to  be  80  per  cent  consisned  at  ejection.  This 
is  close  to  the  predicted  performance;  however  we  should  be  able  to  come 
closer  to  the  firing  specifications  by  slightly  modifying  the  design. 

Let  us  consider  the  general  problem  of  determining  the  required  grain 
design  modifications.  We  shall  make  the  following  assunptions: 

(1)  The  form  function  is  quadratic; 

(2)  Burnt  occurs  near  the  end  of  stroke; 
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(3)  Tbs  loading  density,  A  -  C^lpUgt  ie  '*niall''v  i*  «•»  A  <  0.1; 

(4)  _  is  a  natbeaHitieal  Mulinim. 


AB8iiii>tiGos  (2) ,  (3) ,  and  (4)  axe  met  in  moat  CAD  applications  and  are  not 
felt  to  be  serious  restrictions.  As  for  assuaption  (l)  —  most  fom  functions 
are  cvibiea  but  can  generally  be  rery  closely  approximated  by  quadratic  e:q>res- 
aion.  At  any  rate,  tbe  assuaptions  are  made  for  practical  reasons;  they  lead 
to  rather  slsple  (and  apparently  accurate)  correction  fonnulas. 


We  shall 


y' 


a 


Pm 

c 


hare  need  of  the  following  equations  obtained  from  Section  X: 


=  i'j. 

12lf  ^max 

(£.16) 

(E.  17) 

fJi 

(E.1B) 

w 

=  (^y  Firx, 

(E.19) 

4 

n  +  1  '^max 

(E.20) 

y'*\ 

=  i  i 

f '  \  2  f '  / 

r  -  ^1 

(E.21) 

.  .  ,x  *  X.)  •  J 

(E.22) 

=  Xj  I  +  \,  i’  . 

(E.23) 

Using  these  equations  we  wtuill  derive  egressions  idiicb  deteimine  the  grain 
design  modifications  necessary  to  produce  given  corrections  in  and  v,. 


Differentiating  equation  (E.20)  we  obtain 


(E.24) 
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idiere 


y'  /  1 


(E.25) 


Hie  assiBQitiaa  of  nmali  loading  density  pennits  us  to  treat  the  initial  free 
TOluaa,  as  being  relatirely  independent  of  small  changes  in  the  grain 
design,  i.  e. ,  0^  a  constant.  Hierefore,  equation  (E.22)  maybe  differenti¬ 
ated  to  yield  approximately. 


v'  d  y‘ 


(E.26) 


where 


[i  -  (fj;  -  *  '■'] '  1 

■y'  =  1  . 


{E.'27) 


There  are  essentially  three  grain  design  parameters  that  detemine  ballistic 
perfoimance:  the  web  length,  Lf,,  and  the  two  foim  function  coefficients 
and  Xj  .  We  consider  the  effects  of  bolding  any  two  of  them  constant  and 
varying  the  third. 

Case  1  -  Xj  ,  Xj  constant;  variable 


The  web  length,  L^,  enters  equations  (e.16)  -  (E.23)  only  indirectly 

through  its  effect  on  the  charge  weight,  C{,,  and  the  consequent  effect  on  the 
initial  free  volume,  .  In  the  case  of  small  loading  densities,  we  can  then 
say  that  for  constant  X|  and  X2  a  small  change  in  Lf,  produces  no  appreciable 
change  in  perfoxmance. 


Case  g  —  Xj  ,  Lj,  constant;  Xj  variable 

Differentiating  equations  (E.16)  and  (e.19)  we  obtain 


d  P'y  dp 

=  (3-2.) 


max 


and, 


dy'  =  [y 


{E.28) 


(E.29) 
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CcBA)lnlog  eqiiatlons  (E.a4)»  (E.28)>  axkl  (E.29)  then  leads  to 

f  Jb»  = _ _  . 

^UBX  3  -  2n  +  2  /i  (^  -  f)  (1  -  n) 


(E.30) 


Replacing  differentials  with  finite  differences  the  equation,  (E.30), 
eiqpresses  the  approximate  fractional  change  in  due  to  a  fractional 

change  in  Xj  . 

Similarly,  differentiating  equation  (E.  17)  and  ccnbining  the  result 
with  equations  (E.26),  (E.24),  {E.28),  and  (E.30),  we  obtain 


3-2«  +2/1  {y-y‘)  (1-  n)  Xj 


(E.31) 


which  relates  the  change  in  ejection  velocity  to  a  small  change  in  Xj  « 


Case  3  —  Xj  ,  constant;  X|  variable 

Again  differentiating  equations  (E.  16)  and  (E.19)  we  obtain  in  this 


case 


=  (.-a-) 


max  _  ‘*^1 
c  -r““  » 


max 


(E.32) 


and, 

d  f '  =2  (’y  -  f')  (1  -  n) 

Thsse  equations,  together  with  equation  (E.24)  ,  then  lead  to 

^  ^mex  _ _ 2 _ 

^iMot  3-211 +2 /X  If- f ')(!-*)  x', 


dp 


max 


max 


(B.33) 


(E.34) 


Proceeding  as  in  Case  2,  we  also  obtain 

tJUs.  r  lb  . 

3-  2n-+2/i(f-f')  (l  -  ") 


(E.35) 
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■HiB  results  of  the  three  eases  are  sumarized  in  the  following  table. 

Table  9  iOmniU/aX  CHANGES  IN  UAXHOI  HfflSSORB  AND  EJECTIGN 
YELOCnY  DDE  TO  SMALL  CBANGES  IN  GEAIN  EESIGE  PARAMEESIS 


1  per  cent  increase  in: 

Percentage  increase  in: 

Lhi  (^1  •  ^2 

1^-2  1;*  (^1  »  ^6  coiistant) 

i  (^2  *  constant) 

*Ab8olute  value  of 

^nax 

0 

2 

0 

[*^3-2n)  -m(1-«)] 
j,  |l-2(y-f')(l-«)f2v*-A^ 

where 


r---i  /I 


3  =  1 

A.-*  y 


Ki/z,  r  .=  1 


^'  =  3  -  2ri  +  Zfi  if  -  -y')  (1  -  n) 
r  =  ?-  Z  (^ff  FWK, 

’""1  "  arn~^" 

=  IZ  Ax^/Ui  . 

To  illiistrate  the  use  of  Table  9,  let  us  restate  our  saa^>le  problem. 
Given  a  CAD  and  propellant  with  the  following  properties; 

A  =  4‘.92  in* 

B  =  0.01919  (in*/lbyr)”  in/ sec 
F  =  328,250  ft*lb_^/lb„ 
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n  =0.372 
Uc  =  159. 1  In* 

V  =  12. 185  slugs 
=  13.474  ft 

6  =0  dogrsss 

’y  =  1.295 

p  =  0.05661  lb  in* 

V  =  17.67  in* /Ib^ 

We  wish  to  determine  a  grain  design  to  yield  an  ejeeticn  velocity  of  120  ft/ see. 
Ihe  maxlnun  pressure  must  not  exceed  1600  psi.  Burnt  should  occur  at  or  near 
ejection. 

Bxqploying  the  method  of  this  appendix,  a  tentative  design  is  obtained 
having  the  values 

X|  =  3.1246  lb Id 

Xj  =  17.633  lb„/ in* 

4  =  0.06554  in 
Cf,  =  0.3963  lb„ 

5  =  Xj/p  +  2  Xj  i/p  =  56.20  +  623.0  L  in* 

A  test  firing  of  the  charge  is  made  and  the  following  results  are  obtained: 

=  111.9  ft/ sec 

Pmax  = 

i„,  the  distance  burned  throTi^  the  grain  at  ejection  .=  0.06733  in. 

What  modifications  should  be  made  to  yield  the  specified  performance 
requirements? 
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Solution 

lie  obserre  that  mist  be  ixicreased  by  about  7.2  per  eent  without 
increasing  p  by  more  than  12.8  por  cent.  Using  Table  9>  we  then  proceed 

IHBX 

as  follows. 

Calculate  the  tsIuss: 


-  2  (^-)f  f  Hlj  ■  0.7930 


V',  =  --sr^-<—  =  1.829 

"1  aBfXj  ■ 


12  is. 

T.  =  — TT—  =5.230 


=  y-Vi  [i  *  -/vi ^-4-r]  - 

4^... 

/S'  =  3  -  2n  +  2m  {7  -  y')  (1  -  «)  =  2.820  . 
Table  9  then  has  the  nunerical  form: 


=  0.9156 


Table  10  AHWOTMATE  CHANGES  IN  AND  DDE  TO  A 

ONE  Vm  CEWT  CHANGE  IN  EE3IC«  PARAMETERS  -  SAMPLE  PROBLEM 


1  per  cent  increase  in: 

Percentage  increase  in: 

p 

'^maz 

h 

0 

0 

M 

0.1592 

0.2677 

0.7092 

0.1452 
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In  order  to  detemine  the  corrected  Talues  of  X.|  and  Xj  ne  shall  az4>i- 
trarily  set  as  the  desired  ▼alue  of  the  maxiimn  pressure,  ~  1575 'psi. 

Ha  Bust  then  select  new  values  of  Xj  and  Xj  such  that  the  pTeTio\i8ly 
attained  ralues  ■-  1418  psi  and  .=  ill.  9  ft  /  sec,  are  increased  by  11. 1 
per  cent  and  7.2  per  cent  respectively.  Since  for  email  changes  in  X|  and 
Xj  the  percentage  changes  in  ballistic  perfonnance  are  approximately  addi¬ 
tive  we  may  fora  the  two  equations: 

0.7092  +  0.1592  €,  =  11.1 

0.1452  €j  +  0.2677  €j  .=  7.2 

idkere  e|  and  are,  respectively,  the  percentage  changes  in  X|  and  |  Xj  | 
required  to  produce  the  desired  changes  in  and  v^. 

Solving  for  C|  and  Cj  we  obtain 

~  10.4'  per  cent  (increase  in  X|) 

Cj  .=  21.0  per  cent  (increase  in  |  Xj  |)« 

Ihe  corrected  fora  function  coefficients  are  then 

X,  .=  (1,104)  (3.1246)  =  3.450  lb„/in 

Xj  .=  (1.210)  (17.633)  =  21.336  lb  in*. 

Also,  to  reduce  the  amoimt  of  unbiuned  propellant  at  ejection  we  shall  decrease 
the  web  to  0.0675  in.  According  to  Table  10,  this  will  have  no  appreciable 
effect  on  perforaance.  The  new  chanse  weight  thus  becomes 

Cff  ■-  X|  Lff  *  Xj  Lj*  -  0.3301 

In  sumiary,  'then,  the  reconmended  new  design  parameters  are 
Xj  =  3.450  lb„,/in 
Xj  .=  21,336  Ibp/in* 

=  0.0675  in 
Cf,  -  0.3301  Ib^ 

S  60.94  +  753,8  L  in* 
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Tbese  values  detennine  the  requited  srain  goamtij,  Subatituting  the  valuaa 
into  the  equations  cf  i^ppendiz  A»  leads  to  the  eeleetion  of  a  single-perf 
grain  design  with  an  outer  dieneter  of  0.30  in  and  a  perforation  diameter  of 
0.16  in,  all  surfaces  except  the  perforation  being  inhibited. 
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Figure  2A:  V  Versus  X  and  y  (during  burning) 
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Figure  23:  V  Versus  X  and  ^  (during  burning) 
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Figure  3:  P  Versus  V  and  y  (during  burning) 
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Figure  8.  Values  of  Charge  Weight,  C^j,  Versus  Web 
Size,  L^,,  to  Yield  an  Ejection  Velocity  of  60  ft/sec 
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Figure  12:  p,  v  and  t/Xj,  Versus  x/x^  for  Devices  I,  II,  III 


Figure  13  A: 

Initial  Cross  Section  of 
Seven  Perf  Grain 


Figure  13B; 

Cross  Section  at  Splintering 
of  Seven  Perf  Grain 


Figure  14:  Versus 
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NCMHOJaPRE 

Phlts 

A  -  cross-sectional  area  of  piston 


& 

max 

-  maximum  rate  of  change  of  acceleration 

£'e/sec 

B  - 

burning  rate  coefficient 

feXif})  ” 

C  - 

mass  of  charge  burned. 

Cb  - 

charge  weight 

lb* 

Cy  - 

constant  volunie  specific  hast  of  propellant  gases 

ft-lby 

ib;^‘ 

D  - 

diameter  of  propellant  grain 

in 

d  - 

diameter  of  perforation  of  propellant  grain 

in 

E  - 

1/2  * 

piezometric  efficiency  of  CAD  firing  =  - 

^max^« 

— 

- 

1/2  VvJ 

volixnetric  efficiency  =  — „ — ,7”- 

FpUc 

— 

^2  - 

— 

/(a,  6) 

-  "a  function  of  a  and  b” 

— 

F  - 

impetus  of  propellant  =  SRTy 

ft-lby 

lb* 

- 

resistance  force  due  to  friction 

Iby 

i  - 

acceleration  due  to  gravity 

ft  /  sec  2 

Sc  - 

ccHiversion  factor  .=  3S.  174 

lb  „/ slug 

h  - 

length  of  propellant  grain 

in 

K  - 

spring  constant 

Iby/ft 

- 

friction  constant 

— 

1 


.’rnimmmmsL. 


mHPGBTm  17S8 


-  beat-loss  eoostart 


L-taoriaaD.tal  firing  constant  - 

-  naglsetsd  eo-Toliaa  constant  - 

k  -  aq>irical  corrsetion  factor  taken  to  be  ^  — - 

‘3 

k-  -  ballistic  constants  froai  Section  ttt  rariable 

L  -  distance  burned  throuj^  propeUant  grain  in 

Lff  -  propellant  neb  .=  distance  burned  through  grain  at  "burnt"  in 

Lkmtn  ~  propellant  neb  required  to  attain  "burnt"*  at  end  of  stroke  in 

M  -  anlea  of  gas  per  vnit  naaa  of  propellant 
n  -  burning  rate  eiq>ooent 
k  -  apace  arerage  gas  pressure 
Pg  -  pressure  on  base  ot  piston 


P _  -  maximum  pressure  attained  in  CAD  firing 

onx 


-  corrected  maximwi  pressure 


P  -  dimensionless  pressure  rariable 


=  Hi.  p._Y 

12V  \aBPll 


P*  -  dimensionless  pressure  variable  =  (a ' 


Pjf  -  mathamatioal  maximmi  of  P 

9* 

M 

9 

»» 

M 

•fl 


P^,  -  mexiiniii  value  of  P  attained  in  particular  CAD  firing 


-  maximm  value  of  P*  attained  in  particular  CAD  firing 
ballistic  parameter  = 


qj  -  ballistic  parameter  = 


(f»)  moles 
’  “gram 


Iby/in* 

lb  in* 
lb// in* 

lb //  in* 


2 
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-  ballistic  partmster  .= 
q\  -  ballistic  parsBister  .= 
9*  -  ballistic  parsBieter  .= 

9*1  -  ballistic  parsnetar 


Qi 

Qi 

Q't 

R 


ballistic  pai'SBiBtar 
ballistic  parsnetar 
ballistic  parsoetar 
ballistic  parsnetar 


\t}  m 

0- 

=  <7, /a* 

=  q\/a* 

■  qt/o^ 

=  9i/a 


max 

r  - 

S  - 

5  - 

So  - 
H  - 

t  - 
T  - 
Tv  - 
Ob  - 
Or  - 


univarsal  gas  constant  .=  278S 

-  ballistic  parsnetareif-rc — - ■! 

npy*  -  7y  +  9j 


-  1**3 


burning  rata  = 
propellant  burning  surface 

Ch 

arerage  burning  surface  =  -7- 

P^b 

initial  burning  surface 
final  btiining  surface 
time  measured  from  ignition 
space-average  propellant  gas  tenperature 
constant-volume  flsme  tenperature  of  propellant  gases 
barrel  volune  .=  12  Ax^ 
enpty  chamber  volune 


ft-lb  y  gram 
Ib"rao![e'V 


in/ sec 

in* 

in* 

in* 

in* 

seconds 

•K 

“JT 

in* 

in* 
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Oi  - 

initial  free  chanber  TOliane  .=  Uq  -  0^/ p 

in* 

V  - 

Telocity  of  piston 

ft/ sec 

««  - 

ejection  velocity 

ft/ sec 

V*  - 

corrected  ejction  velocity 

ft/ sec 

F  - 

dimensionless  velocity  variable  =  v 

— 

F'  - 

dimensionless  velocity  variable  ^  F 

— 

- 

value  of  F  at  nhich  P  attains  its  matbenatical  maxinaan 
= 

- 

-  value  of  F'  at  idiich  P'  attains  its  mathematical  maxinun 

=  VY 

— 

^-1 

-  value  of  F  at  coafiletion  of  stroke  if  F^  >  F^j 

— 

‘'l.l 

-  value  of  V*  at  completion  of  stroke  if  F^  >  F,i| 

— 

*^6  - 

value  of  F  at  burnt  - 

— 

H  - 

value  of  F'  at  burnt  .=  ^  Vk 

— 

-  value  of  F  at  completion  of  stroke  if  F^  <  Fj^ 

— 

- 

value  of  F  at  idiich  P  attains  its  greatest  value 

— 

- 

value  of  V  at  vdiich  P'  attains  its  greatest  value 

— 

^.1 

PF* 

-  the  smaller  root  of  the  equation  — 

=  (i-V-^i)  -1 

— 

W 

FF* 

-  the  larger  root  of  the  equation  — 

=  (i  -  ^-r-  •'».)■  - 1 

— 

W-  - 

accelerated  mass 

4 


WJ.  1752 


MIL  HBPOBT 


Wj  -  effective  accelerated  neaa  .=  (l  +  /|)  v 
»  -“l/lc 

X  -  displacement  of  piston  ftxm  initial  position 
-  stress  length 

TOi 

J  -  dimensionless  piston  travel  =  x 


lb. 

slugs 

ft 

ft 


-  value  of  J  at  X  .=  x^ 


a  -  burning  rate  adjustment  factor 


B  -  it-ZJ-il-t 
(Ci,/4)* 


(S) 

in* /lb 


n  -  1 


/S'  -  ballistic  parameter  from  Appendix  E 

HR 


y  -  ratio  of  specific  heats  =  1  + 

^  -  adjusted  ratio  of  specific  beats  =  (>"  l)  '•’I 

A  -  ^i\ 

Y  -  ballistic  parsmeter  =  -y  -  2 

7.t} 

r  -  dimensionless  ballistic  parameter  =  ^  * 


r'  ♦  1 


r*  -  dimensionless  ballistic  paranmter  = , 


loading  density  =  C\jl pU^ 


I  for  y'  =  1 


£  -  fractional  increase  in  a  ballistic  paranoeter 

Vi 

8  -  ballistic  parameter 


7j  -  covolime  factor 


in®  /  lb 


m 


5 


launching  angle 


-  fonn  flmction  coefficients 


linear  fonn  functioD  coefficient  .=  Cj^/ Lj, 
quadratic  fom  function  coefficients 


M' 

V 

V* 


p 

<T 


ballistic  parenieter 
ballistic  parameter  from  Appendix  E 
ballistic  paroneter  .= 

ballistic  parameter  from  Appendix  £ 
ballistic  constants  from  Section  XII 
propellant  density 
ballistic  constant  3 

'P. 


r 


max 


dp 


degrees 

lb„/in 

lb,/ in 

lb,/  in* 


variable 


lb  f/  sec* 


B\(f.p-P^\  dp 
0 

T  -  arbitrary  ballistic  parameter  from  Section  JII 

4>  -  ballistic  parameter  .=  v^*  ^ 

^(y)  -  a  function  of  y* 

ci)(y')  -  a  function  of  y' 


In-iec 

arbitrary 


n  -  piston  diemeter 


in 
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CcRinander,  Operational  Test  and  Evaluation  Force 
U.  S.  Atlantic  Fleet,  U.  S.  Naval  Base 

Norfolk  11,  Virginia  1 

Comnander 

Naval  Air  Test  Center 
Patuxent  River,  Maryland 

Attn:  A.  E.  Adams,  Armament  Test  Division  1 

Commanding  Officer 
Naval  Anminition  Depot 
Crane,  Indiana 

Attn:  Librarian  1 

Comnander 

Naval  Ordnance  Laboratory 
Corona,  California 

Attn:  code  552  1 

Code  561  1 

Ccnmander 

Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland 

Attn:  Library  (HL)  1 

J.  N.  Ayres  (ED)  1 


3 
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praimBUnOW  (Contlmiod) 

Conaandli^  Officer 
Nanral  Ordnance  Plait 
Macon,  Qeorgia 

Attn:  L.  Pertach,  Deai^  Bigineering  DiTiaion  -2 

Catimandlng  Officer 
MBral  Ordnance  Test  Station 
China  Lake,  California 

Attn:  Code  4543  (j.  Sbeiaian)  1 

Code  5663  (W.  P.  Ecxintz)  2 

Code  4072  (R.  Snlth)  1 

Code  4552  1 

Code  4506  1 

Code  5562  1 

Officer  in  Change 

Naral  EiplosiTe  Ordnance  Dispoeal  Technical  Center 
Naval  Propellant  Plant 

Indiai  Head,  Maryland  1 

Comnander 

Naval  Research  Laboratory 
Mechanics  Division 
lashingtcai  '25,  D.  C. 

Attn:  C.  H.  Kingsbury,  Ballistics  Branch  1 

Coananding  Officer 

Naval  Training  Device  Center 

Poi*t  Washington,  New  Toric 

Attn:  D.  D.  Mangieri  1 

Coananding  Officer 

Naval  Onderwater  Ordnance  Station 

Newport,  Rhode  Island 

Attn:  Si/M  ^ 

Oonmonding  Officer 

Nanral  Air  Special  Weapons  Facility 

Alburquerque,  New  Mexico 

Attn:  G.  N.  Bell  1 


4' 


DismiBDricw 


CoBmanding  Officer 
Pleat  inny  Arseoal 
Dover,  New  Jersey 

Attn:  Technical  loformation  Section 
Conmander 

Atwy  Rocket  and  Guided  Ifissile  Agency 
Propulsion  Laboratory 
Research  and  Development  Operations 
Redstone  Arsenal,  Alabama 
Attn:  R.  E.  Betts 

Coomander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 
Attn:  WCLEHR 

SCLSFF-30 

Advanced  Systems  Development  Division 
Cleveland  Bieumatic  Industries,  Inc. 
1301  East  El  Segundo  Boulevard 
El  Segundo,  California 
Attn:  G.  T.  Lampton 

Aerojet-General  Corporation 
11711  Woodruff  Avenue 
Downey,  California 

Attn:  T.  W.  Royer,  Ordnance  Division 

Aerojet-General  Corporation 

P.  0.  Box  1947 

Sacranento,  California 

Attn:  Dept.  56-00  (John  J.  Pabish) 


Aircraft  Amaoents,  Incorporated 
Cockeysville,  liaryland 


Jperiofln  Machine  and  Itemdxy  OeafMny 
7S01  Nbrth  Hetehes  Areeue 
Nilest  lUixiois 

Attn:  D.  L.  Areeon«  Maohanles  Baeeaneh  DiTiBlcn 

keaoar  Research  Rcundation 
W  35th  Street 
Chicago  16«  lUinoia 

Atlantic  Research  Ooiporaticn 
Bdsall  Road  and  Shirley  Higbeay 
ALezandria,  Virginia 
Attn:  G.  V.  Bnersan 
V.  H.  Bolter 

Becknen  and  Ihitley,  Incorporated 
Mieeile  Products  Dirision 
985  K.  San  Carlos  Arenue 
San  Carlos,  California 

Beech  Aircraft  Corporation 
9709  East  Central 
Wichita  1,  Bnisas 

Bell  Aircraft 

Dasi^i  Bflginssring  Airorea  Equipannt 
F.  0.  Box  1 
Buffalo,  New  Toxic 
Attn:  F.  J.  Mnhratti 

Beimite  Poader  Oonpany 
Soledad  Canyon  Road 
Saugus,  California 
Attn:  L.  LoFiego 

Breeze  Corporation,  Incorporated 
700  Liberty  Araoue 
Xlhlcn,  New  Jersey 
Attn:  A.  Lardin 


DiamiBPIKW  (Ooptlnaed) 

Dayton  T.  Broun,  Incorporated 
1305  Stroog  Boad 
Copiague,  Long  laland.  New  Yoric 
Attn:  SiplosiTea  and  Boafc  Backs 

ICr.  G.  R.  Nice 

BJaf 

Britiah  Bil>as8y 
mishington,  D.  C« 

Attn:  D.  G.  T.  Oolebrooke 

Chance  Vought  Aircraft,  Incorporated 
P.  0.  Box  5907 
Dellaa,  Texas 

Attn:  C.  C.  Cox  (Armament  Staff,  Gro\;p  53140) 

R.  D.  Henry 

Cook  Research  Laboratories 

Cook  Technological  Center 

6401  Oakton  Street 

liorton  Grove,  Illinois 

Attn:  U.  A.  Broderick,  Project  Bogineer 

Convair 

A  Division  of  General  Dynamics  Corporation 
Port  Worth,  Texas 

Attn:  K.  G.  Broun,  Division  Research  Librarian 

Douglas  Aircraft  Cospany,  Incorporated 
827  Lapban  Street 
£1  Segundo,  California 
Attn:  R.  G.  McIntyre 
Armanwnt  Gixnp 
Hydraulic  Grovp 
Bogineering  Library 

Denver  Research  Institute 
XBaiversity  of  Denver 
Dhiversity  Paric  Station 
Denver  30,  Colorado 
Attn:  R.  B.  Feagin 
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[CW  (Ooptlmied) 


Edo  Corpoz^ion 
IS  -  10  lUth  Street 

College  Point  5^  New  York  1 

Federal  Laboratories  Ineozporatad 
Saltsbiug,  Pe^ylTanla 

Attn:  R.  B.  Reynolds  2 

Hie  Franklin  Inatitute  of  the  State  of  PennsylTania 
Philadelphia  St  PennaylTania 

Attn:  C.  T.  Darey  1 

Ihlton-Irgon  Corporation 
Box  591 

Dover,  NSw  Jersey  1 

General  Electric  Company 

Burlington,  Vermont  1 

General  Electric  Conpany 
lOS  Plastics  Avenue 

Pittsfield,  liassacfausetts  1 

Grmm  Aircraft  Engineering  Corporation 
Betlpage,  Long  Island,  New  York 

Attn:  S.  Coryell  1 

R.  Howell  1 


Hanley  Industries,  Incorporated 
6143  Bartmer  Street 
St.  Louis  14',  Missouri 

Attn:  T.  G.  Blake,  President  1 

Hercules  Powder  Conpany 
Ballistics  Department 
Kenvil,  New  Jersey 

Attn:  D.  S.  Sinpson  1 
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DiannBDTICW  (continued) 


Holez  Incorporated 
2751  San  JXian  Boad 
P.  0.  Box  148 
BoUiater,  California 

Attn:  E.  W.  Place  1 

Horkey-Moore  Associates 
■24660  South  Crenshaw 

Torrance,  California  1 

Johns  Hopkins  University 
i^lied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  S|prlng,  Maryland 

Attn:  R.  E.  Kemelhor  1 

Kaman  Aircraft  Corporation 
Old  Windsor  Road 

Bloonfield,  Connecticut  1 

Lambert  Bogineering  Conpany 
1100  Macklind  Avenue 
St.  Louis  10,  Missouri 

Attn:  J.  T.  Thorp,  Jr.  1 

Lockheed  Aircraft  Corporation 
P.  0.  Box  551 
Burbank,  California 

Attn:  Staff  Engineering  1 

Lockheed  Aircraft  Corporation 
Technical  Data  Services 
Dept.  72  -  75,  Plant  B-1 
Burbank,  California 

Attn:  C.  C.  Butterfield  1 


Lockheed  Aircraft  Corporation 
Georgia  Division 
Marietta,  Georgia 
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PrangBUHOH  (Oontlnued) 


Lockheed  Missile  end  Space  Dirision 
Santa  Cruz  Test  Base 
Santa  Cruz,  California 
Attn:  L.  T.  Giladett 

Lockheed  Aircraft  Corporation 
P.  0.  Box  904' 

Sunoyrale,  Califoniia 
Attn:  J.  K.  Li^tfoot 
Librarian 

Lockheed  Missile  and  9|)aee  Dirision 

7701  loodley  Aranue 

7an  Nbys,  California 

Attn:  J.  Qurskis,  Jr. ,  Dq>t.  36-30 

The  Martin  Company 
Baltimore  Z,  Maryland 

The  Martin  Company 

Orlando,  Florida 

Attn:  1|/att  Joiner,  Mail  68 

The  Martin  Ccnpany 
Research  Library 
P.  0.  Box  179 
Denrer,  Colorado 
Attn:  Acquisition  Dnit 

McCormick-Selph  Associates 
P.  0.  Box  6 
Hollister  Airport 
Hollister,  California 
Attn:  F.  LaHaye 

McDonnell  Aircraft  Corporation 
P.  0.  Box  516 
St.  Louis  66,  Missouri 
Attn:  Bigineering  Llbz'ary,  D^t;  644- 
V.  W.  Drexelius 


Dlb'IHlBUnClT  (Coatlmed) 

Uiller  Research  Laboratories  DivisiOD 
lliller  Uetal  Products,  Incorporated 
2215  Russell  Street 
Baltimore  30,  Maryland 

Miije  Safety  i^pliances  Cco^any 
201  North  Braddock  Avenue 
Pittsburgh  8,  Pennsylvania 
Attn:  W.  T.  Gofer,  Jr. 

Modem  Metal  Crafts  Coo^any 
222  Diamond  Street 
Philadelphia  22,  Pennsylvania 

National  Aeronautics  and  ^ace  Administration 
1520  8th  Street,  N.  W. 

Washingtcm,  D.  C. 

National  Northern  Corporation 
Vest  Hanover,  Massachusetts 
Attn:  Vasil  Philipcbuk 

North  American  Aviation,  Incorporated 
Columbus  16,  Ohio 
Attn:  Chief  Librarian 

North  American  Aviation,  Incorporated 
9150  E.  Inperial  Highway 
Downey,  California 

Office  of  Technical  Services 
Department  of  Cotmerce 
Washington  25,  D.  C. 

Olin  Matheison  Chemical  Corporation 
Winchester  Western  Division 
East  Alton,  Illinois 
Attn:  W.  B.  Drakes 
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DISmiBUriCW  fOontlnued) 

Olin  Uathieson  Chonical  QorporatioD 
M^  Harai  4v  Oonnectieut 

Attn:  A.  S.  Oogan  1 

Ordnance  Assoc iatea.  Incorporated 
645  El  Centro  Street 

South  Pasadena,  Califoraia  1 

Ordnance  Engineering  Asaeciates,  Incorporated 
407  South  Dearbom  Street 

Chicago  5,  Illinois  1 

Pacific  Scientific  Coapany 
10642  Placentia  Avenue 

Anaheim,  California  1 

Pelmec 

Division  of  Quantic  Industries,  Incorporated 
1010  Connercial  Street 
San  Carlos,  California 

Attn:  I.  W.  nalland,  Ch.  fiigineer  1 

Propellex  Chemical  Corporation 
P.  0.  Box  187 
Edwardsville,  Illinois 

Attn:  E.  H.  Williams  1 

Reacticm  Motors,  Incorporated 
Denville,  N.  J. 

Attn:  Oswald  Williams  1 

Republic  Aviation  Ooiporaticn 
Farmingdale,  New  Yoiic 

Attn:  R.  E.  Fidoten  1 

Republic  Aviation  Corporatism 
223  Jericho  Turnpike 
Hineoia,  New  York 

Attn:  H.  H.  Shea,  Missile  Systems  Division  1 
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MIL  HEPCBT 


DTSTOIEOTIOW  (Oontinued) 


RockBtdyne 

Caeoga  Park,  California 

Attn:  C.  Fingezkood,  Dept.  S91-  368 

Ryan  Aeronautical  Ooapany 
B^gineering  Department 
San  Diego,  California 

Sandia  Cozporation 
Alburquerque,-Neir  liexieo 

Attn:  K.  A.  Saraaon,  Div.  1283^1 

/ 

Sikorsky  Aircraft  Diriaion 
Dhited  Aircraft  Coiporation 
Stratford,  Connecticut 

^pace  Recoreiy  Systeme,  Incosporated 

140  Oregon  Street 

EL  Segundo,  California 

Attn:  J.  A.  VdllB,  Applied  projects 

Stanley  Aviation  Corporation 
2501  Dallas  Street 
Denver  8,  Colorado 
Attn:  J.  B.  Misner 

Rocket  Power/Talco 
P.  0.  Box  >231 
Mesa,  Arizona 
Attn:  B.  P.  Barnes 

Talley  Industries 
P.  0.  Box  876 
Mesa,  Arizona 
Attn:  G.  E.  Hirt 

Temco  Aircraft  Corporation 
P.  0.  Box  6191 
Dallas  22,  Texas 
Attn:  J.  T.  Brewer 


1 


1 


1 


1 


1 


1 


1 


1 


1 
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mamiHUHCB  (ooptlnued) 

Iblokol  Cbeadeal  Cozporation 
Bunter-Bx^stol  DlTislon 
P.  0.  Box  27 
Bristol,  PeonsylTaiia 
Attn:  R.  F.  Morris 

V 

Thiokol  Cbsmical  Oorporation 
Elktoo  Division 
Elkton,  Maryland 
Attn:  M.  D.  Rosenborg 

IMokol  Cbenical  Oorporation 
Hunter-Bristol  Division 
Levittonn,  Pennsylvania 

united  States  Flare  Corporation 
19701  last  Goodvale  Road 
Saiigus,  California 
Attn:  N.  C.  Eckert 

universal  Match  Corporation 
Box  231 

Marlon,  Illinois 

Universal  Match  Corporation 
Annanent  Division 
4407  Co(^  Avenue 
St.  Louie  13,  Missouri 

Halter  Kidde  Conpany 
675  Main  Street 
Belleville,  Row  Jersey 
Attn:  C.  Morgan,  Dev.  Bog. 

Hestinghouse  Electric  Corporation 
Air  Crew  Bogineering  Library 
Mail  Stop  417 
p.  0.  Box  746^ 

Baltimore  3,  Maryland 
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MSimHU'riCW  (continued) 

HBstioftousn  Electric  Coorporatioc 
Adranced  Slystcns  Rigineexdzig 
^ndy  ATenue 
Sunnyrale,  California 
Attn:  Mr.  E.  J.  Barakauskaa 
Ur.  T.  K.  Hoy 

Hydro-Space  Technology  Ineoiporated 
Clinton  Road  and  Route  46 
Heat  Caldmll,  New  Jersey 
Attn:  Library 

Local: 
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